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SUMMARY 

The  Observatorio  San  Calixto  (La  Paz,  Bolivia),  with  the  sponsorship  of  Grant  AFOSR-87-0311  A,  has 
been  studying  Lg  wave  propagation  across  South  America.  Previous  results  may  be  found  in  the  F'nal  Technical 
Report  GL-TR-89-0273. 

Under  Grant  AFOSR-89-0532A  the  earlier  work  has  been  expanded  to  include  other  regional  phases;  in 
particular  the  phases  Li  and  Rg  have  been  investigated  for  a  part  of  South  America,  that  is  to  say,  for 
earthquakes  originated  in  Peru  (or  near  its  coast)  recorded  at  La  Paz  LPB  station.  In  the  future  other  origins 
and  recording  stations  will  be  studied. 

Celedonio  Tito  has  completed  his  thesis  about  Rg  waves  for  Licence  degree  in  the  Universidad  Nadonal  de 
San  Agustin,  at  Arequipa,  Peru.  His  results,  together  with  those  concerning  Lg  and  Li,  appear  in  the  paper 
■Regional  Phases  Lg,  U  and  Rg  across  Peru-Bolivia',  proposed  for  publication  in  the  Bulletin  of  the 
Seismologlcal  Society  of  America,  and  considered  a  part  of  the  present  report. 

Another  kind  of  regional  phase  appears  as  a  coda  of  P  waves,  probably  originated  in  P  wave  scattering, 
possibly  with  S  to  P  conversion.  They  are  remarkable  in  some  Southern  Chile  and  Argentina  earthquakes.  So 
those  phases  are  initially  studied  for  Southern  Chile  and  Southern  Argentina  earthquakes.  The  initial 
considerations  are  introduced  in  the  second  part  of  this  report. 


R.R  Cabre,  E.R.  Minaya,  I.J.  Alcocer,  R.R.  Ayala  (1989)  Propagation  and  Attenuation  of  Lg  Waves  in  South 
America.  GL-TR-89-0273,  ADA218853. 
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The  tectonic  complexity  of  the  Peru-Bolivia  region  is  well  known.  The  fragmentation  of  the  Pacific  plate 
which  gave  rise  to  the  present  Nazca  plate  could  modify  the  subduction  regime,  but  subduction  itself  both 
predates  and  postdates  this  fragmentation. 

The  Arica  elbow  (bend  in  the  coastline  near  18°  -  19°  S)  points  out  a  contortion  of  the  whole  continent 
related  to  the  tectonics  of  the  Andes  and  to  conditions  of  subduction.  At  the  moment  we  do  not  intend  to  analyze 
if  those  relations  imply  a  causality  or  what  kind  of  relations  exist. 

In  this  study  earthquakes  located  in  Peru  and  beneath  the  Pacific  Ocean  and  recorded  at  the  La  Paz  (LPB) 
station  are  examined.  They  have  an  epicentral  distance  range  of  2  to  18  degrees.  Hypocenters  are  restricted  to 
the  following  region: 

3.44-S  through  17.87-S 
68.10°W  through  80.80°W 
Focal  depth  less  that  225  km. 

From  a  preliminary  survey  of  seismic  records  the  zone  of  interest  has  been  divided  into  three  subzones,  in 
agreement  with  local  geology  (see  fig.  1): 

Subzone  I.  Coastline  and  ocean.  !n  the  north  the  trench  is  filled  with  soft  sediments;  in  the  south  a 
crystalline  basement  is  apparent. 

Subzone  II.  Cordilleran  region  (consisting  of  three  geological  provinces:  Western  Cordillera,  Altiplano, 
Eastern  Cordillera).  In  the  west  volcanism  prevails,  except  the  Coastal  batholith  extending  from  6°S  to  16°S; 
the  Eastern  Cordillera  is  made  mostly  of  young  batholiths  and  sediments  folded  and  faulted  (mostly  during  the 
Mesozoic),  but  it  has  not  yet  been  decided  if  that  part  originated  above  an  ancient  subduction  zone  or  was  caused 
by  the  presently  continuing  subduction;  in  the  middle,  south  of  10°S,  a  previous  deep  valley  has  been  filled 
with  detritus  arriving  from  both  Eastern  and  Western  Cordilleras. 

Subzone  III.  Subandean  strip.  It  is  constituted  mostly  of  thrusted  (consequently  folded  and  faulted) 
continental  sediments,  between  the  Cordillera  Blanca  and  the  Brazilian  Shield. 


Seismic  records  obtained  from  the  standard  WWSSN  station  LPB  have  been  analyzed  (see  fig.  2), 
considering  the  possibility  of  a  subsequent  comparison  with  similar  stations  operating  in  other  parts  of  the 
continent.  A  few  records  obtained  from  instruments  of  different  types  located  near  LPB  have  been  compared 
with  those  of  LPB  (see  fig.  3). 

Earthquake  origin  time,  coo,  Jnates,  focal  depth  and  magnitude  have  been  accepted  as  published  by  the 
International  Seismotogical  Center  or  by  U.S.  Geological  Survey. 


Bath  (1956)  identified  the  Li  Phase,  where  i  stands  for  the  intermediate  layer,  in  the  lower  crust;  so  it 
propagates  as  continental  guided  waves. 

Its  apparent  velocity  is  about  3.79  ±  0.07  km/s.  It  is  equivalent  to  the  S*  phase  of  near  earthquakes. 

Payo  Subiza  (1960)  studied  Lg,  Rg  and  Li  Waves  for  earthquakes  recorded  at  the  Toledo  Station.  He 
concluded  that  "the  Li  phase  travels  along  a  continental  path,  especially  through  mountain  path;  the  phase  loses 
energy  quickly  when  traveling  under  ocean  zone,  because  of  the  thin  intermediate  layer.  The  velocity  is  about 
3.84  ±  0.05  km/s. 


Koridalin  (1961)  analyzed  the  U  waves,  finding  a  very  complex  particle  motion.  It  appears  less  developed 
than  that  of  Lg,  but  its  periods  are  longer  than  Lg  periods;  its  velocity  is  about  3.8  ±  0.9  km/s 

According  to  Panza  and  Calcagnile  (1975),  "Li  could  be  used  as  a  sufficient  but  not  necessary  discriminant 
for  structures  containing  a  low  velocity  channel  in  the  upper  mantle*. 

In  70%  of  the  cases  examined  in  this  study  U  is  larger  on  the  horizontal  components  that  on  the  vertical 
one.  But  generally  Li  appears  emergent.  This  prevents  a  precise  time  reading,  so  that  the  velocity  may  not  be 
calculated  precisely.  See  in  fig.  4  the  arrival  time  of  U  vs.  focal  distance.  An  earlier  arrival  generally  means 
a  deeper  focus:  for  a  focal  depth  greater  than  120  km,  individual  apparent  velocity  measurements  vary 
between  3.93  and  3.96  km/s;  for  a  focal  depth  less  than  120  kn  they  are  between  3.79  and  3.86  km/s. 

Li  amplitude  is  apparently  not  influenced  by  depth.  LPB  records  from  offshore  foci  close  to  the  coast 
(distance  less  than  0.5°)  look  the  same  as  records  from  continental  Foci.  However  earthquakes  occurring  on 
the  other  side  of  the  Nazca-South  American  convergent  plate  boundary  do  not  show  any  U  at  LPB.  Probably  the 
contrast  between  oceanic  and  continental  crust  is  a  barrier  for  those  waves.  The  3-component  analog  records 
of  four  earthquakes  were  digitized  at  a  rate  of  3.8  samples  per  second.  These  data  were  used  to  draw  the 
particle  motion,  which  appeared  to  be  transverse  to  the  direction  of  propagation  with  some  horizontal 
polarization  (see  fig.  5).  That  motion  is  clear  for  deeper  earthquakes;  less  so  for  surface  foci,  probably 
because  of  noise  created  by  scattering. 

The  same  digitized  data  were  analyzed  through  the  fast  Fourier  transform.  The  amplitude  of  the 
predominant  frequency  contrasts  strongly  with  that  of  other  frequencies,  both  increasing  and  decreasing  (see 
fig.  9);  this  suggests  that  the  layer  channeling  those  waves  is  more  regular  along  the  entire  zone  intervening 
in  the  transmission  than  would  be  expected  in  that  complex  region. 


Before  the  Li  coda  is  finished,  a  stronger  wave  initiates,  interfering  with  Li,  generally  with  a  clear  but  not 
impulsive  beginning,  often  with  several  arrivals,  but  lasting  a  short  time. 

Press  and  Ewing  (1952)  identified  the  shott-period,  large  amplitude  Lg  phase  as  shear  type  surface  waves 
associated  with  a  continental  structure. 

Bath  (1954)  observed  that  the  Lg  phase  was  composed  of  several  arrivals.  Oliver  and  Bath  (1957,  1958) 
proposed  that  Lg  is  constituted  of  higher  modes  of  Love  waves.  Kovach  and  Anderson  (1964),  Knopoff  et  al 
(1973)  and  Panza  et  al  (1975)  presented  additional  evidence  for  this  interpretation.  Panza  et  al,  (1975) 
demonstrated  that  ‘low  velocity  layers  in  the  crust  are  not  required  for  the  existence  of  the  Lg  wave*. 

The  mean  apparent  velocity  of  Lg  is  3.58  ±  0.012  km/s  (Alcocer  1989)),  which  is  equal  to  the  S  velocity 
in  the  upper  crust  (transverse  polarization  also  favors  consideration  of  Lg  as  S  waves). 

The  amplitude  is  a  function  of  magnitude,  epicentral  distance  and  local  conditions,  for  both  Lg  generation 
and  transmission.  Lg  amplitude  has  been  divided  in  this  study  by  the  corresponding  P  amplitude  to  obtain  some 
measure  independent  of  both  magnitude  and  distance  The  mean  amplitude  normalized  in  this  fashion  is 
3.39  ±  0.426  (a  wide  scatter  of  results  reaching  extreme  values  of  0.2  and  12.0  suggests  some  inadequacy  of 
that  normalization,  but  it  compares  favorably  with  the  Bath  method  related  to  energy). 

For  other  parts  of  South  America  the  normalized  Lg  amplitude  may  be  as  low  as  1 .22  for  paths  from 
Ecuador  to  LPB  and  as  large  as  5.20  for  paths  from  Brazil  to  LPB.  It  has  been  noticed  that  in  most  cases  Lg  is 
clearer  and  larger  at  the  LPB  station  than  at  any  other  station. 
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A  preliminary  Q  value  of  208  for  the  Lg  path  from  Peru  to  Bolivia  was  calculated  from  the  equations: 

A«K*D'1/3  (sin  D)'1/2exp  (-yD) 

Q-  W/(v  yT) 

where: 

A-  amplitude 
D-  epicentral  distance 
y>  anelastic  attenuation  coefficient 
v-group  velocity 
T»  predominant  period 

The  poor  efficiency  value  of  208  (in  the  general  neighborhood  of  200)  appears  characteristic  for  Lg 
propagation  along  the  Andes.  In  contrast,  along  the  Subandean  strip  that  is  covered  with  sediments  Q  is  around 
400,  and  in  the  South  American  Shield  areas  around  900. 


Press  and  Ewing  (1957),  Bath  (1954),  Gutenberg  (1955),  and  Panza  et  al,  (1975)  defined  Rg  as  short 
period  surface  waves  associated  with  continental  structures.  Apparently  it  is  the  fundamental  mode  of 
Rayleigh  waves  in  the  crust,  including  a  low  velocity  surface  layer;  particle  motion  is  that  typical  of  Rayleigh 
waves,  the  group  velocity  being  about  3.0  ±  0.7  km/s. 

Kafka  and  Reiter  (1987)  analyzed  the  Rg  phase  to  determine  characteristics  of  the  Appalachian  crust. 
Tectonic  Characteristics  of  Pg  Phase: 

Subzone  I.  This  seismic  area  is  characterized  by  the  null  or  poor  development  of  the  Rg.  Ocean  earthquakes 
have  no  Rg  apparent  on  the  records.  This  result  is  independent  of  distance  to  LPB  or  focal  depth.  Coastal 
earthquakes  may  have  Rg  poorly  developed,  in  general  emergent  (only  2  clear  Rg  arrivals)  with  period  to 
1.5  s,  with  amplitude  decreasing  for  either  distance  or  depth  increasing. 

Subzone  II.  In  general  Rg  is  clearly  recorded,  independent  of  the  path  of  the  waves;  for  the  northwestern 
Peru  earthquakes  Rg  is  attenuated,  but  until  now  it  has  been  impossible  to  distinguish  if  this  attenuation  is 
related  to  the  epicentral  distance  or  to  any  other  regional  characteristics,  such  as  focal  depth.  The  wave  period 
is  clearly  shorter  than  that  of  the  subzone  I,  that  is  to  say  between  0.8  and  1.2  s. 

Subzone  III.  Rg  arrival  is  always  emergent.  Recording  is  much  clearer  than  for  subzone  II,  (the  same  as 
for  Lg  phases)  independent  of  the  epicentral  distance  or  focal  depth.  Wave  period  is  0.8  to  1 .3  s,  with  a  few 
exceptions. 

Rg  was  recognized  in  121  earthquakes,  generally  with  emergent  beginning,  some  start  times  really 
dubious,  only  a  few  times  with  impulsive  beginning. 

Apparent  velocity  (focal  distance/total  path  time)  varies  from  3.0  to  3.4  km/s.  (See  fig.  7). 

Wave  amplitude  changes  as  a  function  of  distance  and  of  focal  conditions.  The  seismic  waves  arriving  from 
Southern  Peru  have  the  largest  amplitudes.  Those  from  ocean  earthquakes  are  not  visible. 

Records  from  12  characteristic  earthquakes  located  in  different  parts  of  the  region  have  been  digitized;  the 
data  obtained  were  analyzed  in  order  to  draw  particle  motion  and  so  distinguish  wave  type.  It  was  found  that 
particle  motion  is  elliptic  and  retrograde,  a  motion  characteristic  of  Rayleigh  waves,  (see  fig.  8) 
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The  same  digitized  data  were  analyzed  through  the  fast  Fourier  transform;  see  resulting  spectra  in  fig.  9. 

It  may  be  remarked  that  corner  frequency  appears  to  be  between  .8  and  1.06  Hz  for  the  events  in  this  study. 

As  stated  above,  Rg  waves  have  characteristics  of  Rayleigh  type;  they  do  not  appear  to  be  a  simple  mode,  but 
rather  an  interference  of  several  higher  modes,  together  with  some  coda  remnant  of  body  waves.  They  are, 
indeed,  guided  waves. 

Rg  is  practically  non  existent  for  oceanic  earthquakes.  This  observation  supports  the  hypothesis  that  they 
are  transmitted  along  the  granitic  layer:  they  are  of  very  small  amplitude  when  the  path  is  mostly 
cordilleran,  that  is  to  say,  irregular  because  of  the  roots  of  the  mountain. 

Interpretation 

The  Li  waves  are  shear  waves  transmitted  through  the  lower  crust,  from  velocity  and  particle  motion 
observations. 

This  part  of  the  crust  seems  to  be  quite  regular  in  the  Peru-BoNvia  zone  along  which  Li  waves  are 
channeled. 

The  Lg  waves  are  also  shear  waves,  but  transmitted  through  the  upper  crust,  that  is  to  say,  mostly 
through  the  granitic  layer. 

Being  the  upper  crust  a  low  velocity  layer,  energy  is  concentrated  in  it,  therefore  the  amplitude  of  the 
motion  is  larger  that  that  of  Li. 

Q  (that  is  to  say  Qs  for  the  upper  crust)  along  the  path  Peru-Bolivia  is  about  200. 

Rg  has  an  elliptic  retrograde  motion,  predominantly  in  the  azimuthal  plane;  it  propagates  at  the  same 
velocity  as  Rayleigh  waves. 

Rg  waves  interfere  with  the  coda  of  previous  phases;  so  they  are  noisy.  This  makes  any  better  studies 
difficult. 

But  apparently  Rg  is  a  composition  of  several  higher  modes. 
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2.  Map  of  epicenters  and  stations  used  in  this  research.  Country 
borders  are  shown  for  geographical  reference. 
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Part  II 


Research  in  Progress 

In  order  to  study  the  most  relevant  regional  phases  and  how  they  differ  from  classical  phases,  a  systematic 
review  of  South  American  earthquakes  of  magnitude  higher  than  4  that  occurred  during  1978  through  1987 
has  been  undertaken. 

Earthquakes  of  Southern  Argentina  and  Southern  Chile  are  considered: 

37°S  through  52°S 
70°W  through  82°W 

Different  arrivals  within  the  P  coda  were  selected  for  investigation,  after  realizing  that  those  phases  may 
not  be  interpreted  as  Pb,  Pg  or  pP.  Their  study  will  allow  us  to  characterize  earthquakes  and  consequently  to 
relate  features  of  the  seismograms  to  geological  characteristics  of  the  focal  region,  path  and  recording  site. 

A  difficulty  arises  from  the  Influence  of  station  topography  on  the  recording  amplitude,  since  amplitudes 
read  from  the  seismograms  often  do  not  agree  with  those  expected  from  the  epicental  distance  and  magnitude. 
This  difficulty  prevents  an  accurate  calculation  both  of  magnitude  and  attenuation. 

At  this  time  the  diversity  of  residuals  and  ampBtude  within  the  same  zone  is  intriguing,  see  figs.  1 0  and 

11. 


An  east-west  line  at  about  42°S  separates  two  contrasting  zones  according  to  amplitude  and  duration  of 
coda:  records  of  earthquakes  located  north  of  42°S  are  much  longer  than  those  of  earthquakes  located  south  of 
42°S  (see  Fig.  12);  neverless  the  S/Nratfo  appears  similar  for  both  sides  of  parallel  42°S. 

Regional  phases  Lg,  U,  Rg  are  not  apparent  at  LPB  for  earthquakes  with  epicenters  in  this  region. 
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TRAVEI  TIME  (s) 


Figure  10.  Travel  time  curves  of  several  phases  from  Southern 
Chile  and  Southern  Argentina  earthquakes  to  LPB. 
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TRAVEL  TIME  (s) 


500. 


Figure  11.  P  wave  travel  times  from  Southern  Chile  and  Southern  Agrentina 
to  LPB. 
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RESUMEN 


Se  esfcudi An  las  ondas  Pg  originadoz  en  121  sismoc  del  Per 0  <j. 

pr  ift.x  imos  *  su  Costa  registrados  en  La  Pas,  a  la  la::  do  otror: 
estudios  real  in ados  en  otrc-s  par  tcs  del  mundo. 

So  revisan  los  rasgos  general es  do  la  geolog  la  y  la.;  met  odes 
gc.c  f  i  3  icos  apl  icados  al  Area  de  or  i  yen  y  el  r  ec  or  r  ido  de  1  a? 
ondas  Rg;  en  particular  so  consider  an  1  os  per  files  rttr?  pare:'  J< 
estaciones  slsmicas,  anal  is ados  por  James  1971). 

So  han  dividido  los  s i smos  en  estudio  en  tree  son as:  sur 
f.rntro  y  nor  te  deb  ido  a  lor  caractcr 1st  icas  cncentr adas  cn  los 
sismogramas.  La  velocidad  no  t iene  grandes  variants*  una  de  olra 
son  a  entre  2.8  y  3.  A  km/s,  la  pr  imera  zona  present  a  grand. r, 
amplitudes  de  Rg,  sin  diferencia  aprec iable  para  las  diver  sas 
pr  ovine  i  as  geol  Agio  as,  salvo  las  amplitudes  reduc  idas  para  los 
epicentres  en  el  mar. 

En  la  segunda  zona  se  aprccia  una  suave  gradac i6n  en  las 
amplitudes,  quo  s<5n  mayor es  en  la  parte  oriental,  y  disminuyen  en 
las,  provincial  qelAqicas  mas  al  oesfce  sirndo  car  i  inc:<  istentes 
para  si smos  de  origen  oceAnico.  En  la  tcrccv a  zona  las 
amplitudes  son  metier  es  dicha  atenuac i ‘ n  se  debs  al  c ami no 
rocorr  ido,  con  gradac  i«Jn  parecida  a  la  observada  en  el  Peril 
central  (Cordillera  Occidental  y  Cordillera  Oriental'*. 

Una  vez  digit izados  varies  registros  de  diverse  rocorr ido  ’ 
'  "  "  "  -  -i  "-viipi  r-nto  de  par t icula,  teniendo  en 


24 


cucnt  a  la  dcscomposic i An  de  las  ondas  13  (por  ser  el  a-pccto  quc 
mas  purdo  compl  icar  el  anAlisis).  En  cl  mo/imirnto  dt?  part  i cul  a 
prnvaleccn  las  car actor  let icac  de  la  ondas  super f i c i al es  dc 
R-ayl cigh . 

Los  ecpectros  de  Fourier,  obtenidos  con  cl  programa  FTT, 
muestran  una  frecucncia  dc-  corte  entry  D.EJ  y  1.05  !  !_  a  par  t  ir  de 

la  cual  las  amplitudes  disminuyen  r Ap i damente .  Los  mAximos  son 

-5"  -r 

de  0.2  *  10  mm.  a  1.0  *  10  mm.  cor respond iendo  a  frrcuc-nc ias 
de  0.7  a  0.9  Hz.  Los  minimos  inter  cal  ados  no  sAn  dt.  yro 
signi  f  icac  i  An.  Las  ondas  Rg  so  prep  agars  per  la  cortez  a, 
probabl  entente  por  una  cap  a  guia  que  se  encuentra  a  unes  5  A  E  k.;> 
de  profundidad  (capa  granitica)* 
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CAFITULQ  I 


Introduce i6n 


1 .  1 .  General idades 

La  Sismologia  estudia  el  comportamiento  de  la  tierra  como  un 
medio  elAstico  y  tiene  como  final  idad  conocer  ei 
comportamiento  de  las  partes  principaies  de  riuestro  planets 
como  son  la  corteza  ,  man  to  y  nCicleo. 

Mediante  las  ondas  sismicas  que  se  originan  en  el  foco 
de  un  sismo  se  determinan  las  diferentes  caracter 1st icas  del 
interior  de  la  tierra.  Lst-in  clasificadas  en  ondas 
que  atraviesan  las  tres  capas  principaies,  denominadas  ondas 
de  cuerpo,  y  corresponden  a  las  ondas  I'  ,  L,  ,  las  ondas 
super f ic ial es  Love  y  Rayleigh  que  se  propaqan  en  la  capa  mAs 
superficial  (corteza)  y  la  parte  superior  del  manto. 

Nuestro  estudio  radica  en  el  conoc imiento  de  1  as 
caracterlsticaa  y  propiedades  de  las  ondas  super f ic i al es 
Rg,  las  cuales  son  ondas  derivadas  do  la  combination  de  las 
ondas  de  cuerpo  registradas  en  los  sismoqramas  rie  periodo 
corto,  denominadas  ondas  super f ic i al es  guiadas  que  se 
propagan  dentro  de  la  corteza  y  recorren  caminos 
cont inental es  y,  s61o  con  cierta  atenuaciAn,  recorridos 
oceAnicos. 

1.2. -Ob jeto  del  estudio 

El  propAsito  de  este  estudio  es  conoc er  el  comportamiento  de 
las  ondas  Rg  para  di versos  trayectos  del  Peril  a  la  F'az, 
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las  car acter 1st icas  espectrales  de  los  sismos  selecc  ionados 
para  tal  efecto;  como  tambi£n  dar  el  posible  mecanismo  de 
propagac  iAn  de  las  ondas  Rg.  Aportando  de  esta  man  era  un 
tema  de  invest igac i An  poco  conocido. 

Otro  de  los  propAsitos  del  presente  estudio  es  el  de 
presentarlo  como  tesis  a  la  Eacultad  de  Geologia  y  Geofisica 
de  la  Universidad  Nacional  de  San  Agustin  de  Arequipa  ,con 
el  ob.jeto  de  obterier  el  titulo  de  Ingeniero  Geofisico. 

1  •  3.  -Antecedentes 

Press  y  Ewing  (.1352)  t  ident  i  f  icar  on  dos  tipos  de  ondas 
super f ic ial es  de  ba.ja  velocidad  Lg  y  Rg,  ambas  de  recorrido 
continental,  supuestamente  en  la  capa  granitica.  El 
movimiento  determinado  por  Press  y  Ewing  para  la  onda  R9  fur? 
el  ipt  ico  y  retrAgrado  tipico  de  las  ondas  Rayleigh. 

Bath  <1354)  las  define  como  ondas  Rayleigh  de  periodo 
corto  con  una  velocidad  de  grupo  de  3.00  km/s  ,  siendo  6sta 
variable.  Gutenberg  <1355)  da  sugerencias  y  explica  la 
propagaciAn  de  las  ondas  Rg  y  sugiere  que  son  ondas 
atrapadas  en  capas  de  baja  velocidad  en  la  corteza. 

Oliver  y  Ewing  <1357,  1358a,  1358b)  Kovach  y  Anderson  <1964) 
las  interpretan  como  modos  super iores  de  las  ondas  Rayleigh. 
Difcrentes  invest igadores  consideraron  a  la  Rg  como  el 
resultado  de  la  propagaciAn  del  modo  fundamental  de  la 
Rayleigh  a  travAs  de  capas  super f ic ial es,  incluyendo  una 
capa  de  baja  velocidad. 

Anderson  y  Dorman  (1973)  en  base  a  explosiones 
estudiaron  las  ondas  Rayleigh  de  corto  periodo  <0.2  a  1.5  s> 
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sobre  la  geologia  local  de  los  alrededores  de  la  ciudad  de 
New  York.  El  resultado  obtenido  demostrO  la  cxistoncia  de 
una  fuerte  refracciOn  lateral  causada  por  el  contraste  ontrc 
las  rocas  cristalinas  metamOrficas  de  una  parte  de  Manhatan 
con  alta  velocidad  y  la  cuenca  TriAsica  de  Newark  con  ba.ja 
veloc idad. 

Me  Evilly  y  W.  Stauder < 1965)  determinaron  que  existen 
diferencias  entre  las  veloc idades  de  grupo  de  las  ondas  LR 
de  periodo  corto  provenientes  de  sismos  naturales  y  sismos 
ar t i f ic iales,  los  cuales  fueron  registrados  provenientes  dc 
diferentes  recorridos  en  el  Area  del  sol evant amien to  Oz.-rl 
Uplift  y  la  cuenca  de  Illinois,  donde  se  hizo  un  buen 
control  del  sistema  estructural  basamento-sedimentos, 
correspondiendo  veloc idades  de  2.5  km/s  para  pro fund idades 
de  0.65  km,  para  la  caps  semi -in f in  it a  de  3.5  km/s  y  en  la 
parte  de  la  cuenca  del  Illinois  para  una  profundidaJ  de  0.4 
km  le  corresponde  una  velocidad  alrededor  de  1.9  km/s,  para 
0.85  km  de  profundidad  una  velocidad  de  2.8  km/s  y  para  la 
capa  semi -in f in  it  a  que  se  tratA  de  3.5  km/s. 

Alan  Kafka  y  Edmund  Reiter  (1987)  estudiaron  las 
veloc idades  de  grupo  de  las  ondas  Rayleigh  de  periodo 
corto  (Rg)  y  el  rango  de  periodos  estimados  en  este  estudio 
fue  de  0.4  a  1.6  s  . Determinaron  velocidades  de  grupo  de  3.C 
a  3.2  km/s  para  recorridos  paralelos  a  la  gran  estructura  de 
los  Apalaches  y  2.6  a  2.7  km/s  para  recorridos  tr ansver sal es 
a  los  Apalaches,  sugiriendo  que  la  cortesa  deba.jo  de  los 
'  aches  es  lateral  merit  e  anisotropic  a. 
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Capitulo  1 1 


Sintesis  Geol 6a  ica  y.  Geof Isica  del  Area  de  estudio 

1 1 . 1 .  Ub icac i 6n  del  Area  de  estudio 

La  zona  de  estudio  se  halla  ubicada  dentro  las  coordenadas 

geogrAficas  siguientes: 

Latitud  03.442  sur  a  17.87  2  sur 

Long  it ud  68.1  2  oeste  a  00.8  2  oeste 

La  extension  del  Area,  abarca  2f 262, 500  km  de  super ficie. 

1 1 . 2 .  - Qeoloaia  del  Area 

Como  marco  de  referencia  aceptaremos  las  ideas  recopiladas 
por  SuArez  et  al.(1983)  sin  pretender  discut ir las,  aunque 
sin  duda  hay  puntos  dudosos;  pero  es  de  mayor  interns  para 
el  estudio  la  segunda  parte  de  este  capitulo. 

En  la  region  de  estudio  se  pueden  consider ar  cinco 
provincias  geol Agio as  Csegtin  la  apariencia  superficial)  y 
son:  La  LLanura  Costera,  Cordillera  Occidental,  Altiplano, 

Cordillera  Oriental  y  el  Subandino  (fig.  1). 

LLanura  Coster  a 

En  el  norte  del  Peril  la  llanura  costera  es  angosta,  Arida  y 
no  excede  de  40  km  de  ancho.  Limita  al  oeste  con  la  linea 
de  costa  y  al  este  con  la  Cordillera  batol itica  consistente 
pr inc ipal mente  de  rocas  volcAnicas  y  sedimentar ias  de  edad 
Mesozoica. 

En  el  sur  del  PerCi  la  planicie  costera  consiste  en  un 
fuerte  basamento  rocoso  cristalino;  el  PrecAmbrico  costero 
impera  en  la  parte  sur  del  Pen!  llamado  el  bloque  Arequipa, 
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tiene  una  edad  de  1.8  -  2.0  miles  de  mill  ones  de  afios 
(Cobbing  et  al . ,  1377,  Dalmayrac  et  al . ,  1377),  fue  sujeta 
a  defer  mac  i  ‘in  en  la  era  PrecAmbica,  Paleozoic  a  y  Mesozoica. 

Cordillera  Occidental 

Al  oeste  la  Cordillera  estA  const  it  uida  pr  inc  ipalmente  de 
roca  volcAnica  y  plutAnica  del  Mesozoico  y  Cenozoico 
respect ivamente;  tambiAn  se  encuentran  deptisitos  marines  dc 
edad  Mesozoica.  Estas  formas  irregulares  y  continuas  son 
paralelas  a  la  costa.  El  batolito  de  la  costa  se  extiende 
desde  62  a  162  S  paralelo  a  la  costa  (fiq.l);  la  Cordillera 
plut6nica  estA  compuesta  desde  gabros  a  sienita  granitica 
(Cobbing  y  Pitcher,  1972);  t ienen  de  100  a  300  m.a. 

El  plutonismo  parece  ser  joven  y  puede  verse  un 
incremento  en  contenido  silicico;  las  rocas  volcAnicas 
Cenozoicas  se  encuentran  en  la  parte  oeste  de  dicha 
Cordillera,  donde  se  encuentra  una  gran  masa  de  ignimbritas 
que  se  presents  en  el  NeAgeno  (Dalmayrac  et  al . ,  1980). 

En  el  Pen!  central  la  actividad  volcAnica  tiene  su 
inicio  hace  11  mi 11  ones  de  anos  y  par 6  abruptamente  durante 
los  Ciltimos  5  millones  de  aFIos  (Noble  y  Me  Kee,  1377). 
Este  cese  del  vulcanismo  Cuaternario  se  explica  con  el 
inicio  de  la  subducciAn  de  la  placa  de  Nazca  debajo  de  la 
parte  central  del  Peril.  Al  este  y  norte  de  la  Cordillera 
Blanca  las  rocas  desarrol ladas  en  el  Mesozoico  t ienen 
notables  plegamientos.  Los  plegamientos  y  la  compresitin  en 
f’-ta  reqiAn  son  en  general  de  desl  izamiento  nor  flexion: 
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fueron  empu.jados  al  este  sobre  fall  as  inclinadas 
suavemente  al  ocste  en  el  sobreyacente  esquisto  JurAsico 
quo  estA  formando  flexuras  y  pi eqamientos.  El  estilo 
tectAnico  impl ica  una  separaciAn  de  la  cobertura  Nesozoica 
a  partir  de  un  substrate  antique,  con  acor tarn i cute  ex tense 
este-eeste.  En  la  parte  oeste  de  la  zona  se  enc uen t ran 
roc  as  volcAnicas  del  01 igocene  y  el  Miocene,  habiendo 
incon for midad  sobre  rocas  Mesozoicas  deformadas;  la 
de  for  mac  i  An  terming  al  final  del  CretAcico  y  comienzo  del 
Cenozoico  (Dalmayrac,  1978). 


Bloque  Arequipn 

Cobertura  Volcnnica 
Plutones 

Falla  Inverse 

Areas  Falladns 

(1)  LI  anura  Costera 

(2)  Cordillera  Occidental* 

(3)  Al tipi ano  y  Altiplanicie 
(ft)  Cordillera  Oriental 

(r>)  Snbanillno 

(  Da  I  ma v  far  e(  a  I  I ')«()) 


Fig.  1  Principales  Rasgos  Estructurales  de  los  Andes 
Peruanos  (  Audebaud  et  al.,  1973  ) 
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La  zona  plegada  es  truncada  y  limitada  en  el  oeste  por 
plutonismo  granitico  de  la  Cordillera  Blanca  de  hace  3-12 
m.a.  (Stewart,  197-1).  En  el  F'eri'i  central  las  roc  as 
scdimentar ias  del  Hesozoico  estAn  cubiertas  compl etamente 
de  roc as  volcanic as  del  Cenozoico.  No  se  conoce  la  edad 
de  la  de for mac i6n  de  las  rocas  Mesozoic as.  Pec ientemente 
Vicente  et  al.(1979)  mapearon  una  falla  de  ernpuje  con 
direcciifm  al  este  con  un  despl  azamiento  minimo  de  15  I  m 
ccrca  a  Arequipa  (fig.,  1).  En  las  prox imidadcs  a  la  falla 
hay  rocas  antiguas  Pr  ecAmbr  ic  as  sob  re  rocas  sod  intent  ar  i as 
Tlesozoicas;  esta  falla  se  cncucntra  al  oeste  de  la  actividarj 
volcAnica  pero  el  empu.je  ocurrid  durante  el  CretAcico 
tardio,  estando  entonces  al  este  del  arc o  volc&nico  active. 
El  Al  t  ipl ano  y  la.  Al  t  ipl  anic  ie  Central 

La  Cordillera  Oriental  al  este  y  la  Cordillera  Occidental 
al  oeste  limitan  el  Altiplano  y  la  Altiplanicie  del  Peru 
Central.  Al  norte,  alredcdor  de  102  S,  la  Altiplanicie 
Central  desaparece  y  la  Cordillera  Oriental  y  Occidental 
son  adyacentes.  La  anchura  de  la  Altiplanicie  del  Per  i!i 
Central  es  de  10  -  50  km.  de  ancho  y  se  increment a 

consider abl entente  en  el  sur  hasta  casi  200  km.  cerca  riel 
1  ago  Titicaca  en  el  Altiplano.  La  secuencia  va  desde  el 
Paleozoico  a  rocas  marinas  del  Mesozoico  y  se  presenter  en 
el  Altiplano  y  la  Altiplanicie,  aparentemente  se  fueron 
depositando  en  cuencas  prof  Lindas  orientadas  hacia  el  Nor 
Oeste  y  que  despuds  fueron  cubiertas  del  Cretdcico  tardio 
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El  estilo  tectonic*  on  la  Altiplanicie  es  variable.  En  la 
region  de  Huancayo  en  el  Peni  Central,  es  uria  zona  de 
plegamientos  compact  os  y  fall  as  de  empu.je  formando 
cintut  ones  de  10  a  30  kms.  de  ancho,  que  son  separados  per 
zonas  de  suave  pleqamicnto  y  rocas  no  deformadas  de  similar 
anchura  (Leery  y  Davis,  1983). 

En  estas  zonas  de  intensa  deformaciAn  se  pueden 
distinguir  r  ~<cas  del  Paleozoico  o  PrecAmbricas  en  la 
super  fir i"  de  la  estructura.  En  la  parte  sur  del  Altiplano 
estA  cubierta  por  una  gruesa  secuencia  de  molasa 
continental  suavemente  dsformada  depositada  durante  el 
01  igoceno  y  Mioceno.  El  plegamierito  es  menor  que  en  la 
parte  central  y  norte  del  ForAt,  y  las  rocas  del  ferciario 
sedimentario  estAn  generalmente  deformadas  en  plegamientos 
concAntricos  cortados  por  fall  as  inversas. 

A1  nor-oeste  del  lago  Titicaca  las  fallas  de  empu.je 
sugicren  que  la  cobertura  sedimentaria  se  sepani  del 
basamento  subyacentc  y  sc  deformA  durante  el  CrctAcico 
tardio  y  principios  del  Cenozoico. 

Cordil lera  Or iental 

La  Cordillera  Oriental  es  una  zona  extertsa  en  donde 
subyacen  rocas  Pre--liesozoic as  local  izadas  al  ests  de  la 
Altiplanicie  Central.  Hay  masas  cristalinas  PrecAmbricas  y 
rocas  plutAnicas  Paleozoicas  que  afloran  en  Areas  grandes 
del  este  de  la  Cordillera,  par t icul ar mente  en  el  Per  A 
Central.  Las  rocas  sedimentar ias  estAn  const ituidas  por 
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una  gruesd  secciAn  de  estratos  marines  super f ic ial es  y 
cent inontal es  pr inc ipalmente  Paleozoicos  plegados  y 
f r ecuentemento  atravesados  por  fallas  inversas. 

Las  rocas  PrecAmbricas  son  generalmente  rocas 
metased iment ar i as  (filitas  o  esquistos  de  grano  fine).  Se 
encuentran  tambi^n  rocas  cristal inas  del  PrecAmbrico 
(an f ibol itas  granul it icas! .  Se  encuentran  pocas  rocas  del 
Paleozoico  proven ientes  del  Ordovic ico  y  son  general  men  fee 
no  motamir f icas  a  d^bilmente  metamAr f icas.  Las  rocas  del 
PrecAmbrico  superior  y  Paleozoico  son  las  mAs  antiguas 
oxpuestas  on  el  este  de  1  os  Andes  y  son  menos  anisotr  <>p i cas¬ 
tor  mando  una  cortez a  superior  dAbil. 

Suband ino 

La  zona  del  Subandino  estA  const ituida  por  una  franja  de 
rocas  sedimentar ias  plegadas,  paralelas  a  la  Cordillera, 
entre  los  Andes  y  el  Escudo  Drasilero.  El  llmite  entre  la 
Cordillera  Oriental  y  el  Subandino  estA  general  men  t  (a 
formado  por  una  zona  dc  empu.je,  buzando  al  oeste  (fall  a 
inversa).  Las  rocas  estAn  const  ituidas  por  sediinentos 
cont inentales  deposit ados  en  forma  inter mitente  deeds  cl 
Paleozoico  al  Pliocene*.  No  Lay  evidencia  de  un  magmat  i  s.nc* 
And ino.  El  espesor  de  la  cobertura  sediment aria  es  muy 
variado,  no  conoc  i^ndose  6st a  con  detalle.  La  dataciAn  do 
las  rocas  sedimentar ias  del  Tcrciario  Superior  es  pobre  y 
es  dificil  dat ar  con  precision  la  deformaciAn  do  la  zona 

'•’"Hina. 
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La  zona  Subandina  cstA  formada  por  una  series  de 
plegamientos  y  fallas  activas  que  vienen  desde  el  Pliocene 
hasta  el  presente.  La  deformation  se  debe  pr  inc  ipalmente 
al  acortamien to  de  la  cortcza  expresado  por  plegamienlos 
cilindricos  cor t ados  por  fallas  inversas  de  fuerto 
buzamiento  hacia  cl  otiste,  que  penetran  de  las  i  ocas 
sedimentar ias  al  basamento  subyacente  (Megard,1978  ; 
Auderband, 19731 .  Debido  a  la  odad  y  la  intonsidid  de  la 
deformaciOn  tectOnica,  esta  decrees  de  manera  continuada 
hacia  el  este.  Las  sscciones  publicadas  del  Subandino  son 
gencralmente  esquemOt icas,  pues  de  ordinario  no  se  cuenta 
con  dates  de  per  for  ac iones  y  de  per  files  sismicos. 

1 1  ■  3.  -Diver sos  Metodos  Seof  isicos  para  estudiar  1  a  Cor teza  y.  el 
Manto  super ior 

Son  muchos  los  mOtodos  desar rol 1  ados  para  el  estudio  del 
sistema  cor teza  -  manto  superior;  entre  el los  se  pueden 
mencionar  los  siguientess  Metodos  Sismicos,  Qv avimdtr icon, 
y  Magn6ticos.  De  todos  ellos  los  m&todos  sismicos  son  con 
seguridad  los  mOs  directos  y  los  que  muestran  result ados 
menos  ambiguos;  se  usan  los  metodos  de  refracciOn  y 
reflexiOn  de  las  ondas  sismicas  internas  0  de  cuerpo,  para 
el  estudio  del  espesor  y  otras  caracter ist icas  de  las  capas 
que  const  ituyeri  la  cor  teza  y  el  manto  super  ior  .  Entre  otras 
tOcnicas  para  determinar  algunos  rasgos  de  la  corteza 
resalta  el  anOlisis  de  Fourier  con  la  obtenciOn  del 
espectro  de  las  ondas  ilbiett,  1969). 

El  m^todo  de  refracciOn  estO  basado  en  el  cambio  de 
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d ir ecc i An  de  propagaciAn  que  sufreu  las  onda.;  sismicas 
internas  en  el  limite  entre  dos  medios  que  transmiteri  las 
ondas  y  tienen  di farentes  propiedades  fisicas.  Este  mdtodo 
utiliza  las  leyes  de  Snell  y  de  las  velocidades  aparentes 
(Dobrin  y  Olhovich,  1353).  Con  la  aplicaciAn  de  estos 
conoe  imientos,  1  os  sismilogos  pudieron  determinar  la 
estructura  interna  de  la  tiorra  utilizando  1  os  ticmpos 
requer  idos  per  las  ondas  sismicas  para  desplazarse  tlesde  su 
origen  hast a  estaciones  registradoras  situadas  en  punt os 
diet antes,  ondas  que  penetran  en  el  subsuelo,  siguiendo 
trayector ias  que  dependen  de  las  var iac iones  de  la 
velocidad  con  la  profundidad;  este  mttodo  sc.  utilize  en 
todo  el  mundo  para  la  invest  igac  iAn  del  interior  de  la 
tierra. 

El  m£todo  de  reflexion  slsmica  se  basa  en  cl  cheque  de 
la  onda  en  la  super ficie  de  separaciAn  entre  dos  medios, 
donde  parte  de  la  cncrgia  se  refle.ja  regresando  al  medio 
original,  donde  el  Angulo  de  incident:  ia  es  igual  al  Angulo 
de  reflexion.  Este  mAtocJo  css  el  mAs  efectivo  para 
traba.jos  de  prospecciAn  petrolera  pero  tambitAn  se  utiliza 
para  la  invest igac i An  de  la  cortez a  de  la  tierra, 
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Sc?  han  hscho  estudios,  utilizando  el  m^todo  do  Iocs 
espectros  para  deter minar  el  ecpesor  de  la  corteza  y  estA 
basado  en  la  interdependenc ia  que  existe  entre  la 
cstructura  corteza-manto  y  el  compor t am i en t o  espectral  de 
las  ondas  sismicas  que  llegan  a  la  super ficie  de  la  tierra. 
Para  olio  se  calculi  teOricamente  el  efecto  que  producen  en 
el  espectro  de  las  ondas  sismicas  cada  comb inac  i On  de  los 
di fercntes  parAmetros  corteza-manto,  obten iAndose  curvas 
teOricas  y  haciendo  compar ac  iones  con  las  curvas  observadas 
cJe  los  registros  de  una  estaciOn  deter minada.  Este  mfbtodo 
tiene  la  venta.ja  de  determinar  el  espesor  de  la  corteza 
deba.jo  de  la  estaciOn  registradora  en  base  a  r efr  acc  iones  y 
re  flex  iones  segiin  las  velocidades  de  las  ondas  cn  estudio. 

Hannon  (1964),  Phinney  (1964),  FernAndez  S.J.C1965)  y 
Tokuji  (1966),  demostraron  la  apl icab i 1 idad  de  este  mAtodo 
para  determinar  el  espesor  de  la  corteza  y  se  desarrollO  en 
Norte  America;  en  Bolivia  se  realizaron  trabajos  de  Careaga 
(1967)  y  de  Ibiett  (1969).  lambi£n  estos  estudios  de 
la  corteza  y  manto  superior  se  hicieron  mediants  ondas 
super f ic iales  basAndose  inic ialmente  cn  cl  anAl isis  de 
dispersion  de  las  ondas  Rayleigh.  La  dispersion  se  define? 
como  la  variaciOn  de  la  velocidad  de  las  ondas 
super fic iales  con  la  longitud  de  onda.  Gutenberg  (1924), 
fue  el  primero  en  usar  este  mAtodo  en  la  invest igac i On  de 
la  estructura  ter rest  re,  en  particular  las  caracter isticas 
de  la  corteza?  esto  se  logra  en  base  a  la  compar aci On  de 
las  curvas  teOricas  y  observadas  en  los  sismogramas  de  la 
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velocidad  de  fase,  como  tambi^n  para  las  vclocidado:  dc 
grupo.  Los  datos  sismicos  revel  an  que  la  corteza 
continental  tiene  un  espesor  promedio  de  35  km.  y  que  las 
velocidades  promedio  de  las  ondas  longitudinal es  P  y 
tr  diibvei  sales  C  son  de  6.2  y  3.6  km/s  y  sc  van 

incrementando  con  la  profundidad.  En  cuanto  a  su 
compos  ic  i  An  la  cortesa  eat  A  const  ituida  por  tree.  capas 
princi  pales:  la  mAs  exterior  estA  const  ituida  por 

sedimentos  en  gran  parte  consol idados  (capa  sod i ment at  i a ) ; 
su  espesor  promedio  os  de  1  km.,  luego  se  cncuentro  la  capa 
granitica  constituida  por  rocar.  Acidas  con  un  espesor  quo 
fluct&a  entre  10  y  20  km;  finalmente  la  capa  basAltica,  que 
como  su  nombre  lo  indica,  estA  compuesta  de  r  oc  as  bAs ic  as  a 
ul tr abAsicas;  esta  capa  se  prolonga  hasta  la  discont inuidad 
de  liohor ovic  ic .  Los  estudios  de  la  corteza  debajo  do  las 
Andes  los  iniciaron  H. Tat el  y  M.  Tuve  <1956)  y  compart ier on 
un  trabajo  en  la  Carnegie  InstituciAn  de  Washington,  en 
cooper ac i An  con  invest igadores  locales.  Se  hie ier on 
obser vac  iones  do  exploraciAn  en  base  a  exp  1  os i ones  do  minas 
a  una  distancia  de  500  km.;  a  lo  largo  del  fiance  oeste  de 
la  Cordillera  y  Altiplano  no  so  obtuvieron  registros  de  las 
explosiones  ra-As  all  A  de  los  2G3  kms;  esto  se  atribuyA  a  la 
extrema  atenuaciAn  de  la  cortez a. 

Estos  datos  rccolectados  fueron  reexaminados  y 
r einterpretados  por  Wool  lard  (1360)  usando  los  arribos 
secundarios  y  llega  a  un  model o  estr uctur al , que  incluye  una 
4o  la  corteza  do  65  km.  debajo  do  la  Cordillera 
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Occidental.  Estudios  de  Lommitz  <1962)  sobre  dates 
gr avimdtr icos  y  de  explosiones  confirman  el  model o  de 
Wool  lard.  Notemos  que  concluye  que  debajo  del  Altiplano  en 
la  region  tripartita  Chile,  Bolivia  y  Argentina  la  Corteza 
tiene  un  espesor  de  70  km. 

En  1965  una  eerie  de  explosiones  de  las  minas  de 
Toquepala  fueron  observadas  por  estaciones  de  la  red  de 
Arequipa.  El  tiempo  de  origen  de  las  explosiones  y  sus 
coordenadas  son  conocidos  y  se  sabo  que  la  pro  fundi dad 
focal  es  cero.  La  t&cnica  es  descrita  por  James  et  al . , 
<1969)}  los  resultados  fueron  muy  razonables  de  acuerdo  al 
model o  asumido  por  Wool  lard  y  se  admite  una  profundidad 
promedio  de  la  Corteza  debajo  de  la  red  de  estaciones  de  45- 
50  km. 

En  1968  otra  serie  de  explosiones  realizada  por  la 
Carnegie  Institution  se  central izi  en  estudios  de  la  region 
del  sur  del  Per  u  y  nor  te  de  Bolivia.  Los  per  files 
analizados  corresponden  al  sur  y  norte  del  Altiplano  y  la 
profundidad  de  la  corteza  calculada  por  la  curva  de  tiempo 
de  viaje  es  menor  a  60  km.  Estudios  de  explosiones  en  la 
fosa  OceAnica  del  sur  del  Per')  fueron  inter pretados  por 
Fisher  y  Raitt  <1962),  incluyen  un  modelo  de  la  estructura 
de  la  corteza,  indicando  que  la  profundidad  es  alrededor  de 
11  km  debajo  de  la  cuenca  del  Facifico  y  16-17  km  debajo  de 
la  fosa  OceAnica;  este  increment©  de  la  profundidad  lo 
sostienen  deb i do  al  increment©  de  la  velocidad  de  una  capa 
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de  6.7  km/s.  La  corteza  terrestre  de  parte?  de  la  region 
del  Peril  fue  estudiada  por  James  (1371),  utilizando  las 
velocidades  de  fase  y  de  grupo  de  las  ondas  Rayleigh  y  Love 
para  obtener  un  model o  de  la  estructura  de  1  e\  corteza  y  el 
man to  superior.  La  fig.  2  roucstra  1  os  per  files  analizados 
por  James  < 1 37 1 > . 


Fig.  2.  -  Per  f  i  1  cs: 
anal iz ados  por 
James,  1371. 


1 1 . 4.  -Estructura  de  l_a  c or  t ez a  deba.io  del  Per  (  James,  1371) 

Consider amos  el  modelo  de  cortcza  obtenido  por  dispersion 
de  ondas  Rayleigh  de  James  (1371).  Las  trayec tor ia: 
util izadas  se  encuentran  graficadas  en  la  figura  #  2,  este 
modelo  poster iormente  nos  servirA  para  nucstra 
interpretac  iOn.  A  cent  inuac  i  On  daremos  los  model  os  dt? 
cortcza  obtenidos. 

LPB-ARE  <  t  ab  1  a  tf  1  ) 

se-f-l  La  Par-^-r  rtUio  •• 
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ontre  la  mitaci  del  Altiplano  y  mitad  do  la  Cordillera 
CJcc  i dental  . 


TABLA  ft  1 


V 

V 

h 

5.0 

2.8 

2.67 

10 

6.0 

3.  ’16 

2.87 

15 

6. 5 

3.53 

2.97 

45 

3.0 

4 . 45 

3.30 

40 

0.0 

4.40 

n  “• 

O  a  UA. 

140 

3.05 

1. 50 

n  nrr 
lj  •  OvJ 

100 

_ _ 

TCC-CUS  (  tabla  ft  2  ) 

El  perfil 

Toe once- 

Cuzco  (fig. 2) 

cst  A 

reducido  c  as  i 

enter amente 

a  la  Co 

r  d  i  1 1  er  a 

Occ i dental , 

pero  c ubre  una 

corta  secci'in  del 

A1 tipi anc 

,  al 

Sur 

dc* 

Cu^.  c  o . 

TABLA  ft  2 


V 

V 

h 

5.0 

2.88 

2.67 

10 

6.0 

3.45 

2.87 

15 

6.6 

3.  64 

3.02 

40 

7.9 

4.50 

3.  25 

40 

3.  15 

4.62 

n  nc; 
o  a  uJ 

60 

3.  2 

4.65 

3.38 

— 

ARE-NNA  (  tabla  ft  3  > 

El  perfil  Arequipa-waKa  (fig.  2)  se  extiende  a  lo  largo 
do  la  crest a  de  la  Cordillera  Occidental. 


TABLA  ft  3 


V  V  h 


5. 

20 

3.  00 

2.71 

5 

6. 

00 

3.46 

2.87 

20 

6. 

70 

3.69 

3.04 

40 

8. 

00 

4.48 

3.32 

40 

8. 

05 

4.50 

<J  a  u  J 

80 

8. 

10 

4.55 

3.38 

— 
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TCC-ARE  (  tabla  #  4  ) 

Como  el  perfil  entre  Toe once  y  Arequipa  toca  la  plataforma 
continental,  es  probable  que  las  ondas  que  viajan  a  lo 
largo  de  esta  1 inea  scan  refractadas. 


TABLA  tt  1 


Kt 

9 

V 

h 

5.0 

2.  83 

2.67 

5 

6.4 

3.  66 

2.98 

10 

6.7 

3.69 

3.05 

10 

7.9 

4.35 

100 

8.0 

4.38 

o  0(71 
•  uu 

— 

LPB-CUS 

El  perfil  La  Pas -Cuzco  (fig.  2)  se  extiende  a  lo  large- 
del  flanco  oeste  de  la  Cordillera  Oriental. 

El  model o  derivado  (tabla  5)  es  incierto, 
part  icularmente  en  relaci-f-n  a  las  vel oc i dados  del  marito 
super ior . 

TABLA  tt  5 


V 

V 

h 

5*3 

3.05 

2.71 

10 

6.0 

3.46 

2.87 

20 

6.5 

3.58 

2.97 

20 

8.2 

4.68 

3.38 

50 

8.  3 

4.74 

3.40 

50 

8.1 

4.  78 

3.45 

La  fig.  3  es  una  aproximac  i-'n  de  una  secci-in  transversal  a 
lo  largo  de  la  linea  A  A'  (Fig.  2);  cone luyendose: 

1.-  El  modelo  de  tres  dimensi ones  de  la  estructura  de  la 
corteza  muestra  una  var  iaci-f-n  del  espesor  de  11  km. 
(incluyendo  una  capa  de  agua)  en  la  cuencn  oceAnica  a 
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30  km  a  lo  largo  de  la  costa  y  mAs  de  70  km  bajo  la 
Cordillera  Occidental  y  parte  oeste  del  Alptiplano. 


Fig.  3.-  Extructura  de  la  corteza  y  manto  superior, 
James  (1971). 

La  corteza  se  adelgaza  al  Oeste  y  per  debajo  de  la 
Cordillera  Oriental  y  el  espesor  promedio  aproximado 
t’S  de  50  -  55  km. 

La  corteza  bajo  la  cresta  de  los  Andes  se  adelgaza 
hacia  el  norte  y  sur  de  la  region  del  Altiplano,  la 
pro fundi dad  de  la  corteza  maxima  en  esta  parte  de  los 
Andes,  no  asociada  con  el  Altiplano,  <-s  del  orden  de 
55-  60  km. 

2.-  La  vclocidad  de  la  corteza  tiende  a  ser  baja  en  la 
region  de  los  Andes.  Bajo  la  Cordillera  Occidental  y 
Altiplano  hay  alrededor  de  10  km.  de  roca  con 
velocidades  compresionales  y  de  cizalla 

apr ox imadamente  5.0  y  3.0  km/s  respect ivamente. 
Sobreyac iendo  a  una  zona  de  rocas  de  baja  velocidad,se 
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encuent  r a  una  capa  de  ca;i  20  km.  de  profuwdidaJ  dt? 
rocas  con  velocidades  compresionales  y  de  cizalla  de 
aprov;  i  mad  ament  e  6.0  y  3.45  km/5.  Esta  velocidad 

estructural  es  consistente  con  Ian  evidencias  de  la 
cuenca  sedimentaria  del  Altiplano.  La  parte  ba.ja  de  la 
corteza  estA  compuesta  de  rocas  con  velocidades 
compresionales  y  de  cizalla  de  6.6  y  3.6  km/s. 

3.-  Las  velocidades  de  la  onda  de  cizalla  en  el  m:nto 
debajo  de  Areas  cont inentales  tienden  a  ser  del  or  den 
de  4.3  -  4.5  km/ s.  La  mis  ba.ja  velocidad  do  cizalla 

se  determinarA  en  la  sub-cortez a,  debajo  de  la  region 


de  los 

Andes 

en 

1  a 

1  inea  HUA-LPE',  donde 

1  as 

veloc idades 

en 

1  a 

parte  superior  del  man to 

son 

alrededor 

de 

4.: 

28  km/s. 

Esta  region  de  ba.ja  veloc 

idad 

coinc ide 

c 

on 

1  05 

e.je 

s  de  la  alt  a  anomalia 

de 

conduct ividad  elActrica  (Schmucker  ct  al  .  ,  1964-1966). 
Las  velocidades  bajas  se  extienden  a  una  profundidad 
que  no  ex cede  de  los  200  km.  y  el  tope  de  esta  zona 
alcanza  una  profundidad  de  50  -  60  km. 

Ocola  y  Meyer  (1972),  mediante  estudios  de  refraccirtn 
sismica  deter minar <?>n  zonas  do  ba.ja  velocidad  en  la  corteza 
debajo  del  Altiplano  Per Ci-Bol  i via  obteniendo: 

-  La  capa  sedimentaria  metamArfica  de  4-9  km  de  profundidad 
con  velocidades  entre  4.5  -  4.9  km/s. 

-  La  capa  granitica  con  26  -  30  km  de  profundidad  y  6.0 
r-  i  km/s  de  velocidad. 
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-  La  capa  basAltica  de  68  -  713  km  de  profund  idad  y  63  69 

km/s  de  velocidad. 

De  acuerdo  a  los  primer os  arribos  fueron  postal  ados  2  zonas 
do  ba.ja  velocidad  para  la  corteza;  una  entre  los  9  -  12  km 
de  profund idad  y  esta  compuesta  de  mater iales  que  tienen 
velocidades  de  6.0  y  6. 1  km/s,  y  ello  le  atribuye  a  la 
presencia  de  magma  de  cocas  Acidas  volcAnicas  e  intrusivas. 
La  zona  de  baja  velocidad  mas  profunda  sc?  oncuentra 
alrededor  de  los  36  y  46  km  de  profundidad  debajo  del  f'en'i 
y  aprox imadamente  entre  los  30  y  40  km  de  profundidad 
debajo  de  Bolivia,  siendo  atribuidas  a  la  presencia  de 
roc  as  basic  as  volcAnicas  o  intrusivas  y  ello  dc  acuerdo  a 
los  resultados  petrol ogicos  y  geoquimicos  obtenidos  por 
Pichler  y  zeil  (I960),  en  los  ande.:  del  norte  de  Chile. 
Model  os  de  corteza  en  las  tablas  tt  6  y  7. 

TABLA  #  6  TABLA  tt  7 


DEPTH 

VELOC I TY  ! 

<  km) 

(  k  m  /  s 

A 

B 

0.0 

4.50 

4 . 50 

6.4 

4.50 

4.50 

6.5 

6.04 

6.01 

9.0 

6.04 

6.04 

9.  1 

5.00 

5.00 

11.0 

5.00 

5.00 

12.0 

6.  10 

6.  08 

26.0 

6.  10 

6.  12 

28.0 

6.80 

6.75 

31.0 

6.  80 

6.77 

31 . 1 

6.  15 

6.  15 

35.0 

6.  15 

6.  15 

39.0 

6.90 

6.85 

71.0 

6.  90 

6.  90 

72.0 

8.00 

8.00 

90.0 

8.00 

8.00 

DEPTH 
( km) 

VELOC 

(km/ 

I  TY 

s) 

A 

B 

0.0 

4.50 

4.50 

6.4 

4.50 

4.50 

6.5 

6.04 

6.01 

9.0 

6.04 

6.04 

9.  1 

5.00 

5.00 

11.0 

5.00 

5.00 

12.0 

6.  10 

6.08 

30.5 

6.  10 

6.  12 

31.3 

6.  75 

6.75 

36.5 

6.75 

6.78 

36.6 

6.  15 

6.  15 

42.0 

6.  15 

6.  15 

45.5 

6.90 

6.85 

/  0 

6.90 

6.90 

73.0 

8.00 

8.00 

90.0 

8.00 

8.00 
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Capitulo  III 


Car ac  t  er  i  st  i  c  as  de  1  as  Ondas  Rg 

1 1 1 .  1 .  -Descr  ioc  i  <fm  de  1  a  Estac  i 6n  Sismol  Ag  ica  de  LPB 

La  estac i An  sismolAgica  de  la  Pas  -  Bolivia  (LPB)  se 
encuentra  situada  en  la  zona  sur  de  la  ciudad  de  La  Par 
(Ba.jo  Seguencoma)  a  una  altura  de  3,292  m.s.n.m. 

Bus  coordenadas  geogrAficas  son: 

Latitud  sur  162  31'  57.6' ' 

Longitud  oeste  68P  05'  54.1'’ 

La  estaciAn  opera  con  trcs  sismAgrafos  electromagnet i cos 
de  corto  periodo;  cst A  compuesta  de  sismA metros  de  un 
segundo  de  periodo  cenectados  con  gal vanAmst r o  de  0.75  s. 
de  periodo,  la  magn i f icac i An  hasta  octubre  de  3977  de 
estos  instrumentos  fue  de  50,000  veces,  luego  por  o  fee  to:: 
del  ruido  cultural  se  redujo  a  25,000,  con  lo  qut: 
actualmente  se  encuentra  func ionando.  TambiAn  se  tienen  3 
sismAgr  afos  electromagnet  icos  de  largo  periodo  que  conztan 
a  su  vez  de  sismimetros  de  15  s.  y  gal  vanimetros  de 
100  s.  ,  la  magn i f icac i An  con  que  operaban  era  de  1,500 
veces,  reducida  r ec ientemente  a  750  veces.  Esta  estac iAn 
corresponds  al  tipo  WWSSN;  estA  instalada  sobre 
sediment os  Terciarios  de  la  fonnaciin  La  Paz ,  const ituidos 
por  arcillas  consol idadas. 

1 1 1 . 2 .  -  Caracter 1st icas  de  1  as  Ondas  Rg 

Se  han  analizado  121  sismos  entre  los  anos  de  1374  a  1386, 

■  ”  — t  r  ~~  . ■>  '  ■  •  r<f  i  *,n  rini  ,'i  tanto  en 
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el  continents  como  en  cl  oc6ano;  tienen  magnitudes  mayores 
a  4.0  y  azimuts  ep icentr o-estac  i <Sn  entre  los  422  a  1792. 
Los  registros  analizados  pertenecen  a  la  estaei*n  do  LPD  y 
los  dates  han  side  tornados  do  los  bolctines  y  catAlogos 
del  Centro  Inter nac ional  de  Sismolog ia  y  del  U.S. 
Geological  Survey.  La  region  en  estudio  ha  si do  dividida 
en  t-res  zonas,  debido  a  las  car acter  1st  icas  que  presenta 
cada  zona  geogrAfica,  que  podria  mencionarse  a  cada  una, 
como  la  zona  Sur ,  Centro  y  Norte  do  la  region  Peruana, 
adecuAndolo  de  esta  manera  para  un  mejor  reconoc i miento  de 
las  car acter 1st icas  obtenidas  de  los  sismogramas  y  para 
e 1 1 o  menc i onar emos  la: 

1 1 1 . 2 . 1 . -Pr imera  Zona 

Los  sismos  de  esta  zona  (fig.  4)  han  si do  local iz ados 
entre  las  siguientes  coordenadas: 


Lat itud 

Sur 

13.502  a 

18.002 

Long  it ud 

Oeste 

68.102  a 

76.802 

se  analizarfln  35  sismos  (tab  la  If  6)  local  iz  ados  en  toda  la 
region  del  Sur  del  Peril,  con  distancias  epicentrales  de 
1.92  a  8.92  a  la  cst ac  i  '-n  registradora  dc  LPE.<,  con 
magnitudes  entre  los  4.1  y  5.5,  siendo  la  magnitud 
promedio  Mb-4.7.  Tenemos  ademAs  12  sismos  sin  magnitud. 

Las  amplitudes  se  han  consider  ado  en  milimicras  (mu)  y 
las  profundidades  son  variadas  y  van  desde  super f ic ial es 
a  intermedios,  con  un  range  de  variaci<in  desde  los  33  km. 
a  217  km.  A  su  vez  esta  pr imera  zona  se  subdivide  por 
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provincias  geolAgicas,  las  cuales  son:  la  LLam.tr  a  cost  era, 
la  Cordillera  Occidental,  el  Altiplano,  la  Cordillera 
Oriental  y  el  Subandino,  donde  so  registry  un  solo  cismo. 
Sismos  local izados  prAximos  a  la  Costa  (ocAano)  ,  (sismos 
del  1  al  11)  cuyas  ondas  atraviesan  la  LLanurc.  Cost  era,  la 
Cor  d  i  1 1  «=’*'■  a  firr  i  rii^nta)  ,  el  Altiplano  y  la  Cordillera 
Or iental . 

-  Tienen  Rg  con  comienzo  enter  gen  te  (sismos  <■•  ,  0,  10,  11 

salvo  excepc  ional  mento  iisipt.il  etvos  (sismos  1,  7:  fig.  L».  1 .» 

tambien  se  registraron  algunas  fases  dudosas  (sismos  0, 
4,  9)  y  cases  sin  F:g  (sismos  2,  5). 

-  Los  periodor.  predominantes  se  cncucntran  dentro  del 
range  de  1.0  s.  a  1.3  s. 

-  Los  sismos  que  se  encuentran  prAximos  a  la  lines  de 

Costa  (sismos  1,  G,  8,  11)  presentan  un  movimiento  real 

del  suelo  de  740  a  69  mu  a  excepc  i  An  del  sismo  9  cor. 
ba.jo  movimiento  real  del  suelo  de  40  mu;  1  os  sismos  qur 
que  se  encuentran  mAs  al  ocAano  por  ejemplo  1  y  10;  por 
citar  algunos  que  tienen  amplitudes  del  orden  dc:  C3  a 
excepc i An  del  sismo  7  quo  os  mayor  a  los  an ter iores 
mcnc ion ados. 

-  Las  vclocidades  se  encuentran  dentro  del  range  dc  3.11  a 
3.32  km.  /s.  ,  excepc  ional  merit e  un  sismo  con  velocidad  de 
2.85  km/s  (sismo  9). 

--  Sus  profundidades  son  super  f  ic  iales  e  intermed  ias. 

Se  local izaron  cinco  sismos  en  la  Cordillera  Occidental, 
(sismos  del  12  al  16)  los  cuales  cruzan  el  Altiplano  >  la 
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Cordillera  Oriental;  tienen  las  siquientes  carActer ist icas: 
Con  sismos  con  inicio  emergente. 

-  Sus  pro fund idades  son  super  f  iciales  (38  a  91  I  m)  . 

-  Los  per iodos  pr edominantes  se  pncuentran  dentro  del 
range  de  0.9  a  1.1  s. 

--  Tienen  amplitudes  variadas  siendo  los  de  mayor  amplitud 
el  sismo  13  y  el  16: fig.  5.2  <270  y  325  mu)  y  los  de 
baja  amplitud  los  sismos  12  y  14  (80  y  102  mu). 

-  Las  velocidades  se  encucntran  entre  3.28  y  3.01  tun/s. 

-  Se  observan  2  fases  de  Rg  cn  los  sismos  13  y  16. 

-  No  existe  atenuac i un  en  su  recorrido  a  la  estaciAn 
registradora  <excepci6n  del  sismo  15  que  es  atenuado) . 

Sismos  localizados  en  el  Altiplano,  son  los  que  cruzan 

parte  de  la  Cordillera  Oriental  (sismos  del  17  al  30). 

-  Comienzo  emergente  para  los  sismos  18: fig.  5.3,  21,  22, 
23,  y  un  c a so  donde  no  existe  Rg  (sismo  20). 

-  Per iodos  pr edominantes  desde  0.9  a  1.5  s.  y  solamente  en 
el  sismo  27  su  per iodo  baja  a  0.3  s. 

-  Las  amplitudes  se>  encuentran  dentro  del  range  de  60  y 
930  mu . 

-  Las  velocidades  existentes  se  encuentran  entre  los  2.95 
y  3.31  km/s. 

-  No  se  puede  observar  la  Rg  en  el  sismo  20  . 

Se  observan  dos  fases  en  once  sismos. 

-  Las  profundidades  son  super f ic i al es  e  intermedia^. 

Se  ban  anal izado  5  sismos  que  provienen  do  la  Cordillera 

Oriental  y  el  Subandino;  se  caracterizan  por  lo  siguiente: 
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-  Comienzo  emergente  (31,  34,  35.)  y  dec  inpulsivos  (sismo: 
32: fig.  5. 4  y  el  33) . 

-  Periodos  pr edominantes  de  0.9  a  1.2  s. 

-  Las  veloc  idadcs  de  4  sismos  estdin  dent)  o  del  range  cJe 
3.01  a  3.25  km/s.  ,  mas  no  el  sismo  35,  quo  tier..?  una 

veloc idad  de  2.94  km/s. 

-  Tienen  amplitudes  claras  (120  a  930  mu). 

-  Las  pr of und idades  son  super f ic ial es  (sismo  34)  n 
intermedias  para  los  4  sismos  restantes. 

-  Se  observa  s61o  una  fase  de  Rg. 

1 1 1 . 2 . 2 .  -  Segunda  Zona 

La  local  izac  i>in  de  los  sismos  se  oncuentran  dentro  de  1  a~ 
longitudes  72.202  -  80.452  y  latitudes  7.442  13.15  esta 

zona  se  ubica  en  la  parte  central  del  PorQ  (fig. 4),  se  han 
anal isado  63  sismos  qi e  estAn  comprendidos  dentro  del 
rango  de  distancias  de  7.062  a  15.002,  con  magnitudes 
comprondidas  entre  4.0  y  5.6,  para  59  sismos,  exist iendo  4 
sismos  sin  magnitud  <44,  51,  85,  87);  las  pi o fund i Jades 

varian  desde  super f ic ial es  a  inter medios  (ver  tabla  7). 

Se  ha  podido  observar  qua  para  los  si smuts  or  i y i nados 
pr/.ximos  a  la  1  inea  de  Costa  (Oc^ano) ,  corno  los  de  la 
Cordillera  Occidental,  entre  ellos  no  existe  gran 
diferencia,  por  ello  que  los  consider amos  juntos,  los 
cuales  son  11  sismos  (36  al  46)  en  nuestro  anAlisis: 

-  Son  de  comienzo  emergente,  y  los  periodos  predominantes 

pncuentrAn  entre  1.0  1  5 
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•-  Las  amplitudes  son  bajas,  compr end idas  entre  1  os  2B  y 
275  mu.  Be  observA  que  sismos  que  provienen  de  la 
Cordillera  Occidental  no  tienen  Rq  distinguible  (sismos 
44: fig.  6.1,  46),  lo  que  no  sucede  para  lcs  sismos  que 

se  encuentran  pr Aximos  a  la  1 inea  de  Costa  (per  cjcmplo 
los  sismos  40  y  42). 

-  Las  velocidades  estAn  compr end idas  entre  los  3.07  y  3.31 
km/s. 

-  Las  profundidades  son  super f ic ial es  e  intermedia^  de  C  a 
152  km. 

Sr  anal izaron  6  sismos  que  vienen  de  la  Cordillera 

Oriental, los  cuales  tienen  las  siguiontes  car actcr i st ic as  : 
Son  de  comienzo  emergente, 

Los  periodos  pr edominantes  estAn  dentro  del  range  de  1.0 
a  1.2s. 

-  Las  amplitudes  de  los  sismos  varian  entre  20  a  341  mu  ; 

el  movimiento  real  del  suelo  es  mayor  pari  los  sismos 
47,  48:  fig.  6.2  y  49;  siendo  uri  tanto  atenuados  el  50, 

51  y  52. 

-  Las  velocidades  fluctCian  entre  los  2.97  a  3.38  km/s. 

~  Las  profundidades  son  super  f  ic  ial  es  a  intermedias. 

Se  han  local izado  46  sismos  para  la  parte  del  Subandino, 

1 1 cgAndose  a  las  siguientes  cone  1  us i ones: 

-  Comienzo  emergente  y  periodos  entre  0.7  a  1.5  s. 

--  Regular  idad  en  la  mayor  parte  de  las  amplitudes,  muy 
notorias  y  claras  (sismos  72: fig.  6. 3. a  y  93: f ig. 6. 3. b) 
do  22  a  820  milimicras  (mu),  en  algunos  casos  tenues. 
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-  Existen  2  casos  que  no  se  puede  obscrvar  1  a  Rg  '■ '.gsmos 
71  y  97). 

-  En  la  mayor  parte  de  los  sismos  se  pueden  observer 

velocidades  mayores  a  3.0(3  km/s.  y  menores  a  3. 33  km/s, 
con  excepci6n  de  3  sismos  cuyas  velocidades  :c? 

cncuentran  entre  2.93  y  2.99  km/s. 

-  Para  esta  region  se  ha  pod i do  notar  muy  buen  rcgistro  de 

las  ondas  y  podria  ser  deb  1 do  a  que  no  est^n 

in f 1 uenc i adas  per  los  Andes. 

-  Las  profund idades  son  super  f  ic  ial es  e  inter modi  as,  desde 
los  33  a  198  km.  de  profundi  dad. 

1 1 1 . 2 . 3 . -Tercera  zona 

Esta  se  encuentra  dentro  del  rango  de  Latitudes  3.002 
7.302  y  las  longitudes  73.302  -  82.002,  (fig.  4)  en  las 
cuales  se  han  ubicado  los  diferentes  sismos.  Esta  zona 
comprende  la  parte  norte  del  PerCi  y  se  han  anal  uado  23 
sismos  (tabla  8),  comprendidos  dentro  del  tango  de 
distancias  de  9.512  a  18.022,  con  magnitudes  que  estAn 
entre  4.5  y  5.6;  las  pro fund idades  estAn  entre 

super  fic iales  e  intermedias.  Entre  sus  caracter 1st icas 
menc ionaremos: 

Deb i do  a  la  similitud  de  las  caracter 1st icas  entre  los 
sismos  del  oc^ano,  de  la  Costa  y  Cor di Ilea  Occidental  se 
1 es  agrup6  en  uno  solo,  para  el los  diremes  que  son: 

--  De  comienzo  emergente. 

-  Con  periodos  predominantes  variables  dc  1.0  a  1.5  s. 
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7.1.b>  atenuadas  (per  ejemplo  sismo  101)  y  ca:o3  que  no 
se  distinguen  dicha  fase  (cismos  100,  103  y  106). 

-  Las  amplitudes  varian  entre  27  y  302  mu. 

-  Las  velocidades  son  de  3.20  a  3.45  km/s. 

-  Las  profundidades  se  oncuentran  entre  16  a  127  km. 

Para  la  parte  del  Subandino  se  han  analizado  15  sismos 
para  los  cuales  daremos  las  siguientes  car actor  1st icas: 

-  Son  de  comienzo  emergente. 

--  Los  periodos  pr edominantes  son  de  0.3  a  1.1  s. 

-  Las  amplitudes  son  variadas  y  se  encuentran  entre  los 
520  a  24  mu  (sismos  116: fig.  7. 2. a  y  118: fig. 7. 2. b). 

-  Las  velocidades  varian  desde  2.84  a  3.36  km/s. 

--  Las  profundidades  son  super f ic ial es  e  intermedias. 

-  Se  presentan  cases  con  atenuaciAn  de  amplitudes  bajas  de 

24  a  70  mu  (por  ejemplo  sismos  111,  114,  120). 

Es  la  parte  mAs  compleja. 

De  todo  ello  podemos  decir  que  la  Rg  cn  general  son  ondas 
emergentes  pero  existen  algunos  cases  de  tipo  impulsive; 
los  periodos  pr edominant es  estAn  variando  entre  0.8  a  1.5 
s.  las  velocidades  varian  entre  los  2.8  a  3.4  km/s.  las 
mayores  amplitudes  (mi 1 imicr as)  se  han  observado  en  la 

primera  zona  (en  la  1 inea  de  Costa),  las  profundidades 
estAn  entre  super fic iales  e  intermedias. 

En  la  primera  zona  se  ha  leido  la  fase  Rg  con  mucha 

dificultad  para  los  sismos  de  recorrido  oceAnico,  habiendo 

algunas  excepciones  en  que  se  not  a 
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muy  bien.  No  hay 


diferencias  muy  marcadas  en  las  provinc  ias  geolAgicas  quo 
se  han  consider ado. 

En  lar  segunda  zona  se  ha  observado  que  hay  diferencias 
marcadas  en  las  provinc ias  geolAgicas  consider adas  en 
cuanto  a  sus  amplitudes,  siendo  los  sismos  de  mayor 
predominanc ia  los  que  estAn  local izados  en  cl  Subandino. 
Siendo  los  de  mAs  haja  predominanc ia  los  que  vienen  de  la 
1 inea  de  costa  y  la  Cordillera  Occidental. 

La  tercera  zona  se  ha  car actor izado  por  ser  la  de  mayor 

complejidad  en  cuanto  a  la  recopciAn  dc  las  ondas  Rg, 

donde  las  provincias  geologic  as  de  lar  linea  dc  Costa,  la 

Cordillerra  Occidental  y  la  Cordillera  Oriental  son  las 

que  mayor  atenuacitfn  de  ondas  han  mostrado,  siendo  menor 

esta  atenuaci'Jn  en  el  Subandino. 

80"  70° 
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Fig.  4.-  Mapa  de  local izaciAn  de  epicentres  de  las 
tres  zonas  anal izadas. 
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16 
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99 

og 

IV 
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327 
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TABLA  it  8 

(datos  para  la  trrcera  zona) 
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Figs.  5.1  al  5.2.-  Sismogramas  correspondientes  a  la  primera  zona 

Inicio  de  la  fase  Rg. 
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Inicio  de  la  fase 
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Figs.  7.1. a  al  7.1.b  Sistnogramas  correspondientes  a  la  tercera  zona 

Inicio  de  la  fase  Rg. 


Sismogramas  correspondientes  a  la  tercera  zona 
Inicio  de  la  fase  Rg. 


Cap i tul o  IV 


Movimiento  de  particula  de  1  as  ondas  Rq 

IV. 1 . - Descr  jpc i An  del  movimiento  de  particula 

Tara  anal  is ar  c  6  mo  se  muevc  real  men to  el  cubsuolo,  y  on 
particular  de  qu£  tipo  eon  las  onJas  1  .  smic  a? ,  de 
consider an  las  posiciones  succsivas  do  un  pun to  arbitrario 
on  el  interior  do  la  roca,  modianto  la  compos  ic  i  <Jn  do  1  os 
diversos  r eg i sir os  y  su  proyccc  i'in  en  pianos  apropiados. 
Este  proceso  se  llama  anAlisis  del  movimiento  do 
part icula. 

Esta  tdcnica  emplcada  on  la  sismologia,  ha  si do  muy 
b  i en  desarrollada  por  Mutt  1  i  y  IJil  more  <.  l'dGd  .t  y  han 
estudiado  la  relaciin  entre  el  movimiento  do  las  ondas 
incidentes  y  el  observado  en  la  super  ficie  libre.  r-1 
compor  tamiento  do  las  ondas  S  al  llegar  a  una  separ  as  i  >in 
de  dos  medios  diferentes,  por  su  carActer  do  ondas 
t r ansver sal  os  de  torsiin,  results  mAs  complojo  que  ol 
compor tamiento  de  las  ondas  P.  Dado  que  participan  con 
Costas  en  la  qener  aci'in  de  Rgf  consider aremos  de  una  manera 
general  el  reparto  desigual  do  snergia  y  l  a  polar  izac  lit, 
do  las  ondas  result  antes. 

Una  descripci'in  general  del  compor  tarn  ion  to  de  1  as 
ondas  S  vAase  a  continuaciAn:  El  Angulo  do  polarizaci&n 
verdadero  utilizado  on  diferontes  trabajos  do 

invest igac i6n,  est A  dofinido  por  la  ecuac i vn  siguiente: 

Q  -  arc  tg  SH/SV”  < G . 1 > 
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donde  SH  y  S V  non  Ian  comportentes  horizontal  y  vertical 
uel  « w  •-  i  JiZzmicntc  do  1  a  a  ondss  S  ^  pore*  1  ft 

incert idunbre  do  determinar  el  comienzo  de  las  ondas  S,  ya 
f|ue  no  siempre  se  obser  van  claras  e  interfiereri  con  otras 
on  dan,  nos  introducirA  una  cooipl  icac  iuii  en  la 
det or minac  i  An  del  Angulo  aparenteC  ?<  )  en  vcz  del  Angulo 
de  polar  izac  i  An  verdadero  (  ©  ).  El  angulo  de 

polar izac iAn  aparente  dado  por  la  relacifin  siguierste: 

!<  -  arc  tg  SMH/SVH  <  C.2  ) 

donde  SHH  y  SVH  non  las  proyecciorics  en  el  plane- 
horizontal  del  vector  deep las ami onto  de  la  onda  S  <quc  mAs 
adelante  se  detallarA  ampl 1  ament e) .  Para  el  trabajo 
real izado  se  ha  deter minado  el  angulo  de  polar izac i An 
aparente  a  part  ir  de  un  metodo  grAfico.  En  realidad  SHI  I  y 
SVH  son  las  sumas  de  Ice  movimientos  de  la  onda  incidente 
B  y  las  ondas  P  y  S  refle.jadas  en  la  super  ficie. 

En  la  fig.  8  tenemoe  la  representac i An  del  vector 
deep 1 az ami ento  de  la  onda  S  con  respocto  a  un  si  sterna  de 
e.jes  geogrAficos  donde  X  represent  a  el  norte-sur,  Y  la 
direction  esto*  ocr.tc  y  2  el  zenit. 

El  piano  (  XfY  ),  es  el  piano  horizontal  en  el  cual 
proyoctamos  la  direcciAn  de  propagaciAn  del  rayo;  el 
piano  formado  entre  la  direcciAn  de  propagaciAn  del  rayo  R 
(.azimut  del  evento)  y  el  e.je  vertical  Z  es  llamado  el 
piano  de  incidencia  A  piano  longitudinal  vertical.  El 
piano  per pendicul ar  a  la  direcciAn  del  rayo  es  el  piano 
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ijoivlo 


:  iti  tictit ra  el  vector  despl azamiento  dc?  la  onda 


cs  dcnominado  piano  dc  la  onda  C. 
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Si  proyec tamos  cl  vector  despl azamiento  de  la  on da  S 
-obre  el  piano  de  incidencia  y  sobre  el  piano  horizontal 
tendriamos  la  component!?  vortical  SV  y  la  cumponente 
horizontal  SH  r espect ivamente.  Lueyo  si  dure omp onemoc  cl 
vector  SV  proyectando  primer amonte  sob re  la  diiccciAn 
de  propagaciAn  del  rayo  en  cl  piano  horizontal  y  tambiAn 
sob  re  el  e.je  Z  del  zenit  tendr  iamos  las  component""  S  VI I  y 
SVV  del  vector  SV.  Luugo  si  surnames  las  cumponcntes 
SVH  y  cl  SH  obter.emoB  el  vector  resul  tante  SUM  que  es  el 
movimiento  de  particula  proyectado  sobre  el  piano 
hor i z  ont al . 

El  Angulo  de  polar izac i An  verdadero  <@>  estA  formado 
per  el  vector  despl azamiento  do  la  on da  S  y  la  components 
SV  ,  y  cl  Angulo  de  pel  or  izaci  An  aparente  (?<)  cstA 
formado  por  la  component e  SVH  y  la  resultante  SHH.  El 
Angulo  entre  la  direcciAn  del  rayo  incidents  y  la  vertical 
Z  se  denomina  Angulo  de  incidencia  <  i  }  ,  que  viene 
expr  esado  por : 

i  “arc  cost  tg  @/tg  ?<  ) 

El  Angulo  entre  la  direcciAn  de  propagaciAn  del  rayo 
in  ter  sect  ado  con  el  c.je  de  lac  X  ,  nos  darA  cl  Angulo  del 
azimut  del  rayo. 

IV. 2. -Trazado  del  movimiento  de  par t leui a  de  1  as  ondas  Rg 

Para  el  trazado  del  movimiento  de  particula  se  harA  use  de 
las  component es  de  periodo  corto  debidamente  digitizadac. 
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Se  corregiri  c?l  posiblc  desplasamiento  do  la  traza  tomando 
los  valores  de  cada  uno  de  los  registros  mencionados  y 
sacando  el  promedio,  la  recta  que  corrcsponde  a  esc  valor 
se  consider ar A  nivel  cero. 

El  primer  piano  de  los  grAficos  de  movimiento  de 
particula  es  1 1  araado  el  piano  horizontal  y  estA  computes  to 
por  los  valores  de  los  registi  os  N-S  y  E  W,  abaci  sac  , 
ordenadas  v  cspcct i vamentc.  Dicho  grAfico  n,z  Uar A  una 
or i cnt ac i An  prel iminar  del  movimiento  de  particula.  El 
segundo  piano  1  vertical  longitudinal  ')  estA  definido  por  1  * 
dirccc  idn  del  azimut  del  epicentre  (GVH.>  y  la  component  e 
Z  (zenit);  los  valores  de  los  punt os  de  SVH  son  la 
proceed  An  ( su.br  e  la  direcciAn  del  azimut  del  r  ay  o  >  de 
todos  los  puntos  que  ban  dado  or i gen  al  movimiento  de 
particula  del  grAfico  obtenido  en  el  piano  horizontal. 
Sumando  los  valores  de  ambas  componentes  se  obtendrA  el 
segundo  grAfico  corr espond iente  a  ese  piano. 

El  terccr  piano  (vertical  transversal?  estA  competes! o 
por  los  valores  de  la  componente  Z  y  la  proyecciAn  del 
resultado  del  piano  horizontal  en  una  1 inea  perpendicular 
a  la  direcciAn  del  azimut  del  epicentre  que  corresponder A 
al  eje  hor izontal  de  las  abscisas  (SHJ ,  y  mediante  la 
intersecc i An  de  ambos  da  luyar  al  grAfico,  que 
poster i or men te  se  inter pretar  A. 

Rclacionando  los  tres  grAficos,  tenemos  el  movimiento 

•  r*»  ^  v  4-  •!  -  •  1  oc-,  1  -  - 
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podremos  captar  las  caracteristicas  propiar.  dr  las  ondas. 
El  pirograma  para  c-ste  proccso  estA  hecho  eri  el  lenquaje  de 
progr amac i An  basic  y  sc?  1  e  ha  drnominado  qrafmovl  y  ha 
r.ido  ejecutado  per  una  mi cv  omput adore  compatible  (ACER 
-710)  del  Observatorio  Can  Calixto. 

I V -  CJ -  Interpretac  i  An  del  movimiento  de  part  icul  a 

Las  caracteristicas  del  movimiento  de  particula  en  todos 
1  os  cases  most r ados  on  ol  grAfico  del  piano  horizontal, 
ton  un  tanto  irregular  (fig.  9.1. a  al  9.0. a),  pero  se  ha 
observado  una  tendencia  paralela  A  eliptica  a  la  dirccciAn 
del  azimut  est ac  i An-ep icentr o;  a  la  cual  so  le  ha  trazado 
una  linea  arbitraria  (3HH)  que  nos  indica  que  ess 
aprox  imadamente  la  suma  delas  componentes  del  plane- 
horizontal  ello  nos  da  un  indieio  quo  las  particulas  del 
medio  se  est An  moviendo  en  el  piano  vertical  longitudinal. 
Con  la  ayuda  de  este  grAfico  nos  fijaremos  en  los  dos 
grAficos  siguientes;  uno  cuya  proyecciAn  es  paralela  a  la 
direcciAn  del  azimut  y  la  otra  proyocciAn  perpendicular  a 
la  dirccciAn  del  azimut. 

Si  nos  fi.iamos  en  ambe-s,  podemos  observar  que  hay  un  mayor 
predominio  del  movimiento  de  particula  en  cl  piano 
pavalelo  a  la  dirccciAn  del  azimut  (fig.  9.1.L-  al  9.0.b>, 
siendo  en  su  mayor  parte  eliptico  y  retrAgr  ado  lo  cual  nos 
indica  que  estamos  ante  una  fase  de  tipo  Rayleigh  (Bath 
19G1);  En  la  proyecciAn  transversal  a  la  dirccciAn  de 
pr c-payac  i  An,  se  observa  menc-r  predominio  de  components  SH 
(fig.  S.l.c  al  9.3.c).  AdemAs  podemos  observar  que  los 
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Angulos  de  pol ar izac i An  aparentes  dc  loj  registros 
anal izados  para  la  Rg,  se  encuentran  compv  uiididcs  entre 
l,..b  C,  52  '  *  '  -122  .  Csto  no::  set  A  indicando  quo  las  ondas 
son  polar  izadas  do  mode  quo  las  particular  d  c  1  modi.1:*  ce 
estAn  moviendo  en  pianos  vertical  ss  paralelcc  a  1  .a 
dirccc i An  del  movimiento  dc  la  onda,  sin  do.jar  que  ev.ista 
el  movimiento  transversal  a  la  direcciAn  del  <:  n  i  mut . 


70 


as 

CO 


I 


. 

•H| 

O 

u 

fT» 

00 

c 

,3 

OS 

CO 

3 

4-1 

u 

o 

<0 

•H 

Oh 

•H 

m 

00 

XJ 

• 

c 

C 

c 

m 

o 

O 

O 

CD 

H 

*H 

'w' 

<x 

*4 

CO 

*H 

i-4 

✓~v 

0) 

CO 

4-1 

a 

0) 

U 

3 

T3 

•H 

B 

c 

4-1 

•H 

'O 

•H 

CO 

U 

N 

•H 

(X 

H 

a) 

<C 

a 

w 

3 

> 

o 

1 

u 

f-H 

a > 

Cl 

VH 

o 

CO 

CO 

4-1 

c 

o 

4-1 

u 

(0 

CO 

(0 

CO 

rH 

rH 

a. 

w 

a  cu 

0) 

'w' 

rH 

as 

CO 

iH 

TJ 

V4 

CO 

co 

(0 

0 

O 

a 

u 

l—i 

4-1 

*H 

a) 

a> 

c 

c0 

c 

> 

CN 

a) 

\o 

CO 

■H 

c 

•H 

c 

u 

B 

o 

O 

co 

N 

1-1 

N 

o 

u 

< 

> 

•H 

a) 

4-1 

/-N 

O 

M 

4-1 

X 

O 

o 

rH 

3 

x: 

u 

<0 

B 

1 

a 

o 

*H 

• 

o 

V-/ 

*H 

N 

O 

c 

4-1 

< 

ss 

*H 

co 

iH 

J-l 

• 

T— 4 

CO 

cu 

rH 

a> 

Oh 

c 

> 

CO 

o 

4-1 

• 

c 

00  ^ 

CD 

•H 

(0 

> 

Pn 

N-X 

w 

71 


••HO  U\ 
60*0  C  Q) 
pc$  3  CO  (X 

4-»  iH 

CO  *H  Pm  C 
TJ  00  'O 

C  C  /-*  *H 
O  O  O  U 
rH  O 

(0  0) 

rH  H  /-s 

(0  P  o 

QJ  O  o  M 

*o  *h  a  p- 


«0  Pi 
*H  0) 
O  > 
O 

*H  O 
P»  C 
Pi  cQ 
CO  rH 
CL  P-4 

Q)  /-v 

*o  ,£> 


N 

<  rH 
CO 


CO 

4-1 

c 

o 

N 

•H 

M 

O 

.C 

O 

C 

C0 


ON  PX 


> 

o 

22 


CO 
Pi 
0) 

> 

0) 

c 

to  >■"* 

Pi  Pi 

a 

CL  rH  *H 
CO  N 

c  o  < 

'O  -H 
■HUH 

own 
o  <u 
<u  >  u 
>>  to 

O  i — <  i — t 

h  to  3 
CL  O  O 
-H  -H 
■M  -o 
»~t  v-<  c 
to  a)  ai 
e  >  cl 


00  «— N 

■H  to 

Ml  ' — 


72 


Evento  N°  16  Az*  281°  (Estac-Epi)  &**19 


■  -HI  O 

I 

00  *o 

c 

0 

P 

CO 

P 

CO 

4J 

rH 

O 

CSI 

CO  *H 

PH 

*H 

u 

•o  00 

TJ 

m3 

C  C 

/-N 

C 

o  o 

a 

0) 

w 

ex 

/"V 

«0 

»H 

•H 

rH  rH 

4^N 

<0 

ex 

CO 

4J 

ex 

w 

a>  a 

P 

1 

T3  *H 

a 

c 

u 

•H  'O 

cO 

CO  U 

N 

4J 

rH  Q) 

< 

u 

(0 

P  > 

u 

w 

O 

rH 

01 

VI  o 

CO 

*->  c 

o 

M  CO 

cO 

u 

(0  H 

«H 

ex 

H<  Ph 

0) 

rH 

<u  ^ 

TJ  x> 


c  M 
at  4-1 
•h  c 
s  o 

•H  N 
>  t4 
O  M 
X  O 
x. 

i 

.  o 
c 

co  at 

•  *—t 
Ot  fU 


60  - 
*W 
Ex  ' 


M  rH  O 
mm  <r 

MB)  csi 

cx  M  cn 

at 

C  >  H 
to  M  /— v 

•H  c  4-»  N 
a  m  3  C 
u  U  B 
at  4J  -i-i 

NO 
O  rH  <  1-- 
M  M 

a.  o  m 

-H  .-H  O 
4J  Z 

rH  M  M 

M  at  m 

C  >  rH  O 

4J 

§ 

£ 


73 


Fig.  9.4  Sismo  N°  72 
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Evento  N°  78  A z  =  322°  (Estac-Epi) 
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Cap i tul o  V 


An  A 1  i  s  i  ^  He  peipp;  t  r  os  de  I  as  on  das  Rq 

Introduce  i An 

La  apl  icac  i  >in  de  lac  series  de  four  icr  r,r  ha 

intensi f icado  en  el  campo  de  la  geofisica,  como 

c onsec uenc i a ,  cn  parte,  de  las  posibi 1 idades  ofrcc idas 
por  1  o 3  ordena d c r  e s  e  1  e c  t  r  <!•  n  1  c o r».  La  si  s m o  1  o q  i.  a  f  o r> 
pr obab  1  ement o  dentro  de  la  geofisica  la  que  mAr.  interns 
most r 4  deeds  el  principle  per  estas  tAcnicas  y  hoy  c:: 
rasi  imposible  imaqinar  tanto  en  prospecc  i  An  omo  en  la 
sismoloqia  de  terremotos  un  tr  at  ami  onto  de  lor.  dat  os.  que 
no  se  real  i.:t’  part iendo  de  la  apl  icac.  i  An  de  la 

tr ansformada  de  Fourier. 

Esto  obedecc  on  parte  a I  propio  car  Ac ter  de  les 

dates  sisaicos  y  al  tip...  de  pr  oblema  e.  .■  el  lor.  prerentan . 


fintre  1 

0  5  tpilia: 

i  inAs  trecuontc.nente  t  r  at  ados 

c L  A  n 

los 

est  ud i os 

do 

dispersion  y  atenuaciAn 

do 

on  d  a ' 

super f  ic 

ialer.  y 

el  pr oblema  inverse  de  la 

estruc fcura 

c  or  t  r:a- 

(Dcin  t  o  f 

las  i nver t i gac i ones  d  e  mo 

■:  x  5fti'. '  / 

■.  aractcv  is 

t  i 

car  focal  os,  la 

ostructura 

del 

n  dc 1  to , 

cl 

estud  io 

de 

la  respuesta 

de  los 

instrument's 

de 

r  eg  i  st  r  o, 

et 

c.  MA abundante 

es  a  tin  la 

apl  i 

caciAn  de 

1  05 

espec  t  r  os 

en  prospecc i An  si  sin  lea  donde 

se 

emp 1 ean 

1  OS 

pr i nc ip ios 

generales  de  la 

transfer  ,iac 

i  An 

de  Tour 

ler  , 

t a mb i An  en  lo  relacionado  con  toda  la  teoria  del  filtrado 
temporal  y  er.pacial,  dcsarrollada  para  el  tratamiento  de 
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serial -ruido 


1  os  dates  con  el  fin  do  mejorar  el  cociente 
basado  en  die  ha  tv ansformac i6n. 

Por  :u.  parte,  en  prospecci4n  gravim£trica  se  ban 
incorporado  t&cnicas  de  anAlieis  arminico  y  fil trade  a  la 
interpretac  i  On  de  mapas  y  per  files  (los  m£todos  de 
prolongaciOn  anal  it  ica  >  seyur.da  der  ivada  eon 
aplicaciones  rle  filtros  especiales).  Tambittn  sc  apl  ice 
en  pr  ospecc  i  On  magnet  ica  para  filtrar  los  mapas 
magnet  ices,  en  el  m6t odo  CE7  en  pr ospecc i On  electrics 
para  determinar  cl  contrasle  dc  las  di ferentes  capae 
subyacentec  y  en  muchos  otvos  c amp os  come  la  oceanogr a f i a 
etc.  Por  todo  cste  trabajo  que  realize  la  tr  e.nsfor  coda  Jc 
Fourier  es  de  particular  importancia  aplicarlo  al  cstudio 
dc  las  ondas  Py,  ye  que  nee  moot  r  ar  0  las  car  ac  lev  i  st  3  •- ...  = 
espectrales  y  el  compor tamiento  de  la  onda  a  tra.tis  do 
los  Andes  del  Per 0 . 

Descr  ipc  i  On  teOr  ica  do  1  a  tr  ansformada  de  F our  i er 
Com  Cm  monte  llamada  la  integral  de  Fourier,  es  introducida 
general mente  c  prrtir  del  desarvollo  en  serin  del  mismo 
nornbre  para  funciones  per  i  Od  icac ,  al  hacer  infinite  el 
intorvalo  que  define  el  peri odo  fundamental  de  la 
funciOn,  obteniendo  luego  la  transformada  de  Fourier  que 
actCta  directamente  sobre  la  tunc  i'!m  original  del  tiempo 
f(t>  para  pasar  al  dominio  de  la  frecucncia. 

Es  necesario  para  el  anAlisis  dc  los  eopect.  00  Je 
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la  integral  do  Fourier,  para  lo  cual  se  darAn  con  raucho 
detail  e  los  principles  bAsicos  dr  la  teoria  rspoctral 
para  luayo  aplicarlo  al  cAlculo  numdr ico;  para  olio  nos 
apoyaveiTios  en  las  series  de  !  uur  ier,  par  a  lucgo  llegar  a 
la  transformada  do  Fourier. 

Si  consider  amos  una  funciAn  f(t)  per  i  A  cl  i  c  a  do  u.ri 
per  iodo  (T>  quo  puede  represent  arse  per  la  ::or  io 
t; r  i gonomtj 1 1  i c a  en  td-rmi nos  clc  suma t  or  i  a : 

f(t>  ~  a«/2  +  ta^'-os  nw,t  +  b^sen  nw,t)  (7.1) 

_ <*  =  L _ 

tomando : 


m 

tq  z  -  -  ;  w  „  =  2  TT/  7 

—  '  O 


la  cual  medianto  procesos  do  sustituciAn  de  variables  y 
pr oced imientos  tr igonomfitr icos  obtonemos: 

f(t:»  -  c„  y~Yi  c^co-<n  <7*2> 

'ft  -  i  _ _ 

llamada  la  eerie  tr igonomdtr ica  do  Fourier,  la  cual  es 
urra  funciAn  peri  Ad  ica,  y  represent  a  la  suma  dc 
componentes  sinusoidal es  que  tienen  diferentes 
frecuencias.  La  componente  sinusoidal  de  frecuencia 
w0  =  2  7T  r  -  2  TT  /  T,  se  1  a  conoce  como  la  frecuencia 

angular  fundamental,  por  tener  el  mismo  per iodo  de  la 
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funciAn  y  ~c:r  cl  primer  anniuico  conocido.  La  componcnte 
w  _  =  n  w_  se  dcnomina  la  e?n£>sima  armAnica  de  la  funciAn 

periAdica  porque  es  en  base  a  esta  variable  que  van 
cambiAndose  las  f recuencias.  Los  coef icientos  c  m  y  1  os 
Angulos  zm  se  conocen  como  amplitudes  armAnicas  y 
Angulos  de  fase  respect ivamente.  • 


Son  import  antes  la:,  propiedades  de  las  funciones 
seno  y  coseno  para  evaluar  los  coef icientos  de  Four  it;  ; 
estas  propiedades  for  man  un  con. junto  de  funcioi.2 
ortogonalcs  cri  el  interval  o  —  T/2  <  t  <  T/2  .  LJtilizando 
la  expresiAn  (7.1)  y  1 ae  funciones  ortogonalcs,  sc 


obtienen  los  coeficientes  de  Fourier  a  p  ,  a  m  ,  b^  : 


T/2 

:/T  \  fCt!’  dt  ; 

-T/2 
T/2 

b _  -  2/T  1  f(t)  sen  (n  w  t) 

"I  I 

-T/2 


Los  c o'  •  f  icicntee  de  Four  i  nr 


funciAn.  Para  hacer  uso  de 
discreta,  es  com  An  exprtisar 
los  exponencialcs  complejos 
per i Ad i c as  fit). 

LI  seno  y  el  cosciio  e 


2/T 


c  os  <  n 


L  )  dt 


dt  ;  para  n  -  1,2 . 

determinan  un 1  vac  ament r  la 

los  espcctros  de  frccuencia 

estas  series  en  t Armines  de 
♦  jnw  t 

a  ,  i  :.rj  funciones 


pueden  cxpr'oar  en  t&rminos 


de  exponencialcs  complejos  : 
Asi  (7.1)  se  convert iria  en : 


f  Ct 


an/2 


+ 


jnw  t 

l/2(am  -  j  b  m  )  c 


1/2  (a  r.jb 


-jnw  t 
)  c 
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que  vendria  a  ser  la  formula  t;7.3>. 
f un c i 6 n  ha r e m o s  q u e : 


->b« 

-  (a  +  jb  > / 2 


Par  a  i  mp  1  i  f  i  c ,  1  a 


Entoncess 


f  ( t  :> 


/n  z  i 


(c- 


■jnwe  t 


utiliz&ndo  las  pr opiedades  dt  las 


-jnw, 1 

i-  c  e  1 

—m 

sumatorias  final  monte 


t  enemor : 


la  cual  se  dcnomina  la  eerie  compleja  cJc  Fourier, 
donde  los  coef ic iontes  c m  y  c  ^  ,  podemoa  cvaluarlos 
f  inalmente  on  t6rminos  de  los  coef ic ientes  de  Fourier  a ^ 
y  b  quo  se  han  dado  con  anterior idad  tenicndo  en 
consider ac i *n  que: 

-  jnw,  t 

e  -■  cos  nw.t  -  j  sen  nw0t 


T/2  -jnw, t 

1/T\  f(t)  e  dt 

-T/2 


6  t  amb  i  d?n : 


c 


<n 


1/T 


f  <t> 

0 


-jnw,  t 
e  dt 


c 

-  o 


=  1/T 


jnwe  t 
e  dt 


n  *=  1,: 


donde: 
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(donde  s  indica  que  cs 
conjugado  c omp  1  e  jo  ’ 


entoncGs: 

~  1/2 

para  valores  n  i*  0  ,  para  lo  cual  c  es  real  y  es: 


■:0  =  ae/2 

la  gvAf ica  del  valor  de  los  coeficientes  ooiiiplc-jos  c 
versus  w  ( frecuenc ia  rngular)  es  dcnomin a  espcctro  de 
amplitud  de  la  func idn  per  iAdica  f  <t ) .  La  qr  if  ica  dr] 
Angulo  de  fase  <  r.  *  de  c  versus  w,  s_  denomina  espcctro 
de  fase  de  f(t>.  Sabemos  qu:.  u  tiene  va  lores  enter  os,  los 
espectros  de  amplitud  y  de  fase  no  son  curvas  continu-c, 
sino  que  aparecen  en  la  variable  discreta  nw;  por 
consiguiente,  ss  les  denomina  espectros  >.!  >  frecuenc  ia 
discreta  o  espectros  de  1  incus.  La  repci  s  :  c.nt  ac  i  An  de 

los  coeficientes  complejos  c  versus  1^  variable 

nwt  , espec i f ic a  la  funciAn  periodica  f(t>  en  el  dominie- 
de  la  frecuenc ia,  asi  como  f(t>  versus  t  especifica  la 


funciAn  en  el  dominio  del  l  i  tempo. 


Las  series  de  Fourier  son  un  instrumento  mu> 
poderoso  en  el  tratamiento  de  di versos  problemas  que 
implicaii  funciones  per i Adicas.  En  el  campo  de  l  a 
invest igac i An  de  funciones  no  periAdicas  ,  es 
necesario  desarrollar  un  mil  tod  o  de  anAlisis  de 
Fourier  adaptado  y  para  ello  sc  introduce  la  integral 
de  Fourier  y  los  espectros  continues. 
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Si  tensmor,  una  func  i  An  per  iidica  f  <  t  >  de  per  iodo 
T  y  so  hace  que  tienda  al  infinite,  cntonces  la 
fund  An  resultants  f(t?~  1  im  fT<  t )  deja  de  ser 
periAdica  .  El  problema  reside  en  encuntrar  la 
rcpresentac i An  de  Fourier  de  csta  func i An  no 
periAdica  y  para  el  lo  tenemos  la  forma  exp  oi.enc  ial  dc 


la  eerie  de  Fourier. 

»<*>  -  Z 

*  »  at 

sabemos  que: 


jnw0t 


T/2  -jnw„  t 
f  <  t  >  e  d  t 

-T/2 


dondc-  ftt)  tambiAn  ~s  puede  oxpresar: 


haccmcc  que  el  T  tienda  a  infinite  cntonces  we 

tiende  a  cero.  Sea  wc  -  A  w,  cntonces  la  frr.cucncia 
de  cualquier  armAnico  nw  debe  cor  responder  a  la 
variable  general  de  frecuencia  que  describe  cl 
espectro  continue  .  Es  decir  n  tiende  al  infinite  a 
medida  que  w„  =  ^w-»0,  obteni^ndose  un  valor  finite, 
ecto  es: 
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n  w0  =  n  A  w  -*•  w 


de  ecta  manera  (7.5 )  ce  convierte  en: 


f  <  t)  - 


f  T/2  jn  a.;,  t  jir4«t 

f(t>  e  dt  e  zjw 

-T/2 


en  el  11 mite  T  —+  o&  ,  A  w  dv  y  la  '..Lima  tori  a  ce 
convierte  en  la  integral  rob  r  a  w,  c:  deiir,  la 
func  i6n  no  periodica  f  <t  >  re  convierte  ms 


<*>  r  ,  «> 


f(.fc)  "  1/27T 


-jwt 

f  (  t  ’  n  d  t 


jwt 

(3 

_ 


J  W  i.  ,  . 


ii  re  define: 


-  jwt 

F  (  W )  ~  f  (  t  >  £?  tJ  t 


veremos  que  (7.8)  se  convierte  en  : 


fCt)  -•  1/21 T  F<w>  dw 


(  7  .  O  > 


La  fun: i An  F  ( w )  def  inida  en  (7.7)  re  cu.occ  como  la 
integral  de  Tourier  o  tr anrformada  invrrna  de  Fourier. 
A  eetar  don  funcionen  m  lan  llama  tambidn  par  de 
tr  nnsformadar.  de  Fourier  .  La  condition  sufic  iente 
pero  no  necesaria  para  que  exista  F  (w>  estA  dado  por: 


:  f  c  t :» !  dt 


teniendo  que  ser  finito  el  valor  absolute  de  la  integral . 
La  funci<f>n  F(w)  es  compleja  y  se  puede  expresar  como: 


F(.w>  -  f  (t  >  (cos  wt  -  j  sen  wt  >  dt 
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F  <w> 


f  <t ) 


f  ( t  >  -en  vt  ri  t 


F  <  w  )  - 
donde  : 


wt  dt 


j  l:  (  w ) 

A(w>  -  j  EKw)  =  ! F (  w  > ! e 


A(w>  ~  parte  real 
i  LKw)—  parte  imaginaria 

:f<w>:  ~  acw>  +  j b ( w :• 

FI  par  de  Iransformadas  tienen  medioc  rquivalcntc:; 
dr>  r  epr  csent  ad  An  :  uno  en  el  dominio  del  tiempo  , 


f  <t > 

• 

f 

y  el  otro  en 

el  dominio 

de  1  a  ft  ecu  end  a, 

F  (  w) 

• 

La  funciAn  (7.7) 

transforma 

la  funciAn  f(t>  , 

que 

se 

encuentra  en 

el  dominio 

del  tiempo,  a  su 

funciAn  equivalents  F(w)  ,  en  el  dominio  de  la 
frecuencia  y  la  ocuac  i An  (7.8!)  invierte  el  process. 

La  fund  An  (7.7)  ,  anal iza  la  funciAn  del  tiempo 
en  un  espectro  de  frecuencia  y  la  ccuaciin  (7.8), 
sintetiza  el  espectro  de  frecuencia  para  obtencr 
nuevamente  la  funciAn  en  t£rir.inos  del  tiempo.  A  1  £< 
fund  An  F(w>,  en  valor  absolute  !f(w>!,  .;c  1  a  denomina 

espectro  de  amplitud  y  I:(w’»  espectro  de  fasc  cJc  f(t)  y  si 
elevamos  al  cuadrado  el  espectro  de  amplitudes  tendriamos 
el  espectro  de  potencia. 


'.’.3.-  Aol  lead  An  del  AnAl  isis  Espectr  al  a  Dndas  Rg 

La  evaluaciAn  numArica  de  la  transformada  de  fourier  a  las 
ondas  en  estudio  que  provienen  del  Perti, 
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se  hizo  mediante 


un  programs  llamado  FFT  (Fast  Fourier  Transform),  hecho  err 
lengua.je  de  pvogramaciin  Fortran  77,  el  cual  fue  a  dec  uado 
para  la  microcomput  adora  Acer  710  del  CJbser  vator  i  o  San 
Calixto;  dichc  proyrama  sc  apcya  en  la  tr  an:;  for  mad  -• 
directs  da  Fourier  y  :;■!.!:■  mediantc  una  ser  ic  de 
procedimientos  art  i  f  ic iosos,  se  logra  increment or  la 
velocidad  del  proccso  y  por  lo  tanto  disminuir  el  tiempo 
empleado,  per  ello  se  la  denomina  transformada  rApida  de 
Fourier;  todo  el  process  del  pi  ograma  est A  represent ado  an 
el  diagrama  de  bloques  de  la  fig.  10. 

La  microcomput adora  nos  per  mite  calcular  ingentes 
cantidades  de  datos  y  petv  a  estos  resultados  se  ha 
utilizado  la  informaci-in  digital;  al  process  de  convert  ir 
los  datos  analogic os  en  digitales  se  le  llama 
digit  isac  i<f*n.  La  funciAn  slsmica  muestreada  a  interval  os 
regulares  no  coincide  exactamente  con  la  traza  inicial;  es 
pues  necesario  consider ar  cl  efecto  que  esta  sustituciAn 
introduce  sobre  los  resultados  del  anil  iris  dc 
f r ecuenc i as;  de  aqu i  que  para  definir  una  f unci  An  mediantc 
valores  muestreados,  es  precise*  tenor  un  minimo  ric  dos 
valores  por  ciclo,  es  dec  ir  el  intervale*  de  muestreo  ha  ric? 
ser  igual  o  menor  que  T/2  siendo  T  el  per  lode*. 

Do  aqui  que,  dada  una  funciAn  fit),  el  espectro  de  la 
funciAn  que  se  obtiene  muestr cAndol a  a  intervales  T  silo 
tiene  valor  para  frecuencias  iguales  A  infer  lores  a: 

w  =  1/C2T) 
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Pig.  10  Diagrama  tie  bloqucs  del  programa  F  F  T. 
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la  cual  se  denomina  frecucncia  dr?  Miquist.  Para  nuc~tro 
estudio  se  consider A  un  inlcs  .alo  de  muestreo  de  0.25  s. 
do  periodo  y  el  proceso  de  d  ig  i  t  izac  i  An  se  real  izi 
mediante  una  digit izadora  model o  6/D.D.  RS-3b  do  la  Gerber 
Scientific  Instruments  Company.  Con ven i en t oment c  so  Man 
digit izado  un  total  de  12  registros  de  periodo  corto, 
cor respond ientes  a  12  sismos  quo  estAn  local irados  en  las 
tres  zonas  del  Area  de  estudio  con  el  propAsifco  de  obtonor 
un  mayor  conocimiento  de  las  car  aeter  ist.  icas  cspec  t  r  a3  cv;  y 
las  pr opiedades  de  las  ondas  Rg,  para  el lo  se  menc ionar  An 
1  os  diferentos  pasos  quo  siqu,.  el  programs  de  espectros: 
CorrecciAn  de  los  datos  digitizados.  Los  dates 
digit izados  se  rncuentran  on  unidadcs  do  digit  izac  i  >!n;  ha 
sido  neccsario  introducir  un  factor  do  corrocciAn  para  quo 
los  mismos  scan  consider ados  on  unidades  mAtricaF  on  el 
cAlculo  de  amplitudes;  ol  valor  del  factor  de  corrccciin 
consider ado  e;  2. 12G5  mm. ,  siendo  esta  la  parte  primera 
que  realiza  dicho  programa. 

EliminaciAn  de  la  deriva  lineal.-  En  la  operaciAn  de 
obtener  la  funciAn  muestreada  es  muy  coimSn  quo  la  linen 
base,  o  1 inoa  cero,  vonga  afectada  por  una  deriva  mis  o 
menos  acentuada;  para  e  .minar  esta  deriva  so  considers  en 
ol  programa  como  linoa  cero,  o  base,  la  definida  por  las 
medias  di  itmAt  icas  de  las  amplitudes,  donde  1  funciAn  del 
tiompo  f  <  t )  ost  A  definida: 

A(I)  =  Amplitudes  de  f<t)  medidas  a  interval  os  fijos. 
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DT  -  Intervalo  de  muestreo. 

Filtro  de  pasa  baja 

Far  a  eliminar  el  ruido  de  frecuencias  alt  as  se  utiliza  un 
suavizado  de  5  puntos,  de  acuerdo  a  las  siguientes 
f  drniul  as: 

« 

-  Para  1  os  dos  primer os  puntos  se  utilize: 

A  <  1 :»  =  (fi(l)+A(2))/2 

AC2)  =  <Aa)+2*A(2)  +  A<3>)/4 

-  Para  1  os  dos  ultimos  puntos  se  utilizes 

A(N  •  1 )  =•-  <A<N-2)+2*A(N-l)+A<N> >/4 

ACN)  =  (ACN-D+ACN)  >/2 

-  Para  los  puntos  intermedios  se  consider^  que: 

existe  un  bucle  que  varia  desde  1=3  hasta  NN=N-2  siendo 
N  el  ndmero  total  de  puntos  a  filtrar,  y  que  se 
almacenarA  en  A(I).  Para  ello  se  utilize  la  sig. dents 
formula:  x 

A  < I )  =  C A  < 1-2) +2*A <1 - 1 ) +3*A  C I ) +2*A (I  +  1 ) +A  < I +2) ) /9 
es  un  filtro  simdtrico  que  supone  que  no  so  produce 
desfasaje  y  es  conocido  como  suavizado  en  pequeFlo, 
dcnominac  i  >in  utilizada  como  contraposic  i>in  al  tipo  de 
suavizado  en  que  los  datos  son  sustituidos  por  un 
pol inomio  de  grado  mAs  o  menos  elevado,  ajustado  por 
minimos  cuadrados. 

CAlculo  de  la  tr ansf ormada 

Despu^s  de  efectuar  el  filtrado  de  pasa  ba.ja,  se  procede 
al  cAlculo  de  la  transformada  de  Fourier,  que  viene  dada 
por  la  exprositfin  siguiente: 
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r  - ° 

Iw)  ~  fCt)  [cos  wt  +  .j  sen  wtj  dt 

J  _  «? 


para  luego  reducirla  a: 


FCw)  = 


A  C  w )  +  .  j  BCw) 


donde: 


AC  w)  = 


fCt)  cos  wt  dt 


BCw)  = 


fCt)  sen  wt  dt 


El  cAlculo  num^rico  de  FCw)  se  realist  con  la  subrut ina 
FFT,  en  la  cual  so  ha  utilizado  el  algor itmo  de  Cooley  y 
Tukey  <  1365)  y  es  una  adaptacitin  de  la  prograniaciin 
desarrollada  por  Clearbout  et  al .  en  196G. 


En  el  programa  la  parte  real  corresponds  a  ACw)  y  la 
parte  parte  imaginaria  a  BCw),  las  cuales  son  calculadas 
por  separado  y  es  a  partir  dc  el  las  que  so  calcula  el 
espectro  de  las  amplitudes,  que  os  igual  a: 

FCw)  =  /aCW)2  +  B  C  W  f 

donde  a  dichos  espectros  so  les  hizo  la  correspondicntc 
correcciin  por  la  curva  de  respuesta  del  instrumento. 
VII.4.-Interpretaci6n 

La  ventana  de  tiempo  utilizada  para  anal izar  la  serial 
sismica-fue  de  32  segundos,  tornado  a  partir  del  tiempo  en 
que  comienza  la  fase  Rg  incluyendo  la  maxima  amplitud  de 
las  ondas.  Par  a  convert  ir  la  serial  sismica  do  la  forma 
anal  igica  a  digital  se  hizo  uso  dc  la  d ig i t  izac  i <* n;  donde 
la  serial  sismica  os  convert  ida  cn  un  n timer o  deter  mi  nado  de 


punt  os  digit  izadoi 


Para  obtener  nuestros  espectros  sc- 
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hi  so  Ltso  de  128  punt  os  de  digit ijaci An,  que  fueron 
obtenidos  do  la  component©  vertical  de  los  si  smog ramas. 

El  intervale  de  digit izaciAn  utilizado  fue  de  0,25  s, 
teniendo  que  aceptarse  un  error  instrumental  do  la 
digitizadora  de  0,0166  s. 

4 

El  range  de  frecuencias  dentro  del  cual  sc  traba.ja  para 
obtener  los  espectros  de  amplitudes  fuA  deterrainado  per  el 
interval©  de  periodos  leidos  en  los  sismogramas;  dicho 
range  de  frecuencias  es: 

0.666  hertz  de  frecuencia  minima  que  os  correspondientc  al 
perl ode  mAximo  de  1.5  s,  leido  en  los  sismogramas  y  1.428 
hertz  de  frecuencia  mAxima  que  le  corresponde  a  0.7  s  de 
periodo  min imo  leido  en  los  sismogramas. 

Se  han  aceptado  algunos  criterios  para  consider ar  la  serial 
sismica  apropiada  para  obtener  los  espectros.  En  los  cases 
en  que  existen  grAficos  de  movimiento  de  particula,  se  les 
hizo  un  seguimiento  y  se  consider  A  la  serial  a  partir  de  los 
puntos  digit izados  que  mostraron  un  movimiento  eliptico 
r etr  Agrado. 

AdemAs  en  todos  los  cases  se  consider*  lo  siguiente: 

-  Para  la  exact itud  del  anAlisis  so  tomaron  una  cantidad  de 

puntos  digit  izados  de  la  serial  sismica  que  sea 

n 

exactamente  igual  al  exponcncial  de  2  . 

-  La  vclocidad  correspondiente  a  la  serial  del  primer  punto 

de  digit izaci An  tornado  para  la  obtenciAn  de  los 

espectros,  estA  entre  los  3.19  a  3.06  km/s  para  los  12 
espectros  anal izados. 
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-  La  velocidad  correspondiente  a  la  sefial  del  (iltimo  punto 
de  digit  i^ac  iiin  no  debe  exceder  de  los  2.8  km/s. 

Para  nuestra  interpretac iin  consider aremos  4  sismos  cuyo 
paso  ha  sido  por  debajo  de  los  Andes  y  8  sismos  cuyo 
recorrido  lo  hagan  por  el  Subandino  y4  Cordillera  Oriental, 
los  cuales  estAn  ubicados  en  la  figura  11. 

De  los  4  sismos  que  atraviezan  los  Andes  se  determinaron 
las  siguientes  car acter 1st icas  espectrales: 

-  En  la  figura  12.1  (sismo  7),  la  mAxima  amplitud  espectral 

fue  de  0.4214  E-5  mm.  y  le  corresponde  una  frecuencia 

predominante  de  0.6875  hertz;  la  minima  amplitud  do 
0.1163  E-6  mm.  y  le  corresponde  una  frecuencia  de  1.1563 
hertz . 

La  frecuencia  de  corte  correspondiente  es  de  1.033  hertz; 
el  espectro  de  amplitudes  nos  inuestra  que  a  partir  do  la 
frecuencia  de  corte,  al  incr ementarse  las  frecuencias,  el 
espectr.o  disminuye  mostrando  una  fuerte  caida  de 
amplitudes  (caida  de  energia). 

-  La  figura  12.2  (sismo  12),  tienc  una  mAxima  amplitud  do 
0.2037  E-5  mm.  y  una  minima  amplitud  de  0.3171  E-6  mm. 
y  les  corresponde  las  frecuencias  pr edorninantec  do  0.750 
y  1.2188  hertz  respect ivamente. 

Su  frecuencia  de  corte  es  de  0.773  hertz;  el  espectro 
nos  muestra  una  fuerte  caida  de  amplitudes  y  por  ende  una 
fuerte  caida  de  energia  a  partir  de  la  frecuencia  do 
corte. 
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-  El  e spectra  de  la  figura  12.3  (sismo  IE!)  tiene  una 

amplitud  maxima  de  0.9714  E— 5  mm.  ,  siendo  r.u  frecuencia 
predominante  de  0.688  hertz,  la  minima  amplitud  es  de 
0.4335  E-6  mm.  y  su  frecuencia  de  0.969  hertz. 

La  frecuencia  de  carte  es  0.906  hertz,  sienda  la  caida  de 
energia  a  partir  de  el  la  similaV  a  los  antcriorcs 
espectros. 

-  En  la  figura  12.4  (sismo  22)  la  amplitud  maxima  fu6  de 

0.8045  E-5  mm.  y  su  frecuencia  predominante  de  0.898 

hertz,  siendo  la  minima  amplitud  de  0.1457  E-6  mm.  y  le 

corresponde  una  frecuencia  de  1.365  hertz. 

La  frecuencia  de  corte  es  0.93  hertz,  mostrando  fuerte 
caida  de  las  amplitudes,  siendo  tambidin  similar  a  los 
anteriores  espectros. 

Para  la  segunda  parte  del  anAlisis  tratamos  8  sismos  fig. 10 
de  mayor  recorrido  del  Subandino  en  relaciAn  con  la 
Cordillera  Oriental,  y  sus  car acter 1st icas  espectrales  lo 
daremos  a  cont inuac i An: 

-  La  figura  12.5  (sismo  93)  nos  muestra  una  amplitud  del 

espectro  maxima  de  0.1844  E-4  mm.  y  su  frecuencia  de 

0.742  hertz,  la  amplitud  minima  registrada  es  de  0. 1852 
E-6  mm.  que  tiene  una  frecuencia  de  1 , 343B  hertz. 

La  frecuencia  de  corte  fue  de  0.915  hertz. 

-  En  la  figura  12.6  (sismo  78)  la  amplitud  mAxima  fue  dc 

0.5275  E-5  mm.  y  su  frecuencia  de  0,719  hertz,  la 

amplitud  minima  calculada  fue  de  0,1643  E-6  mm.  y  su 
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frecuencia  corrcspondiente  de  0.906  hertz. 

Se  obtuvo  una  frecuencia  de  cortc  de  0.93  hertz. 

La  figura  12.7  Csismo  72)  muestra  una  amplitud  espectral 
maxima  de  0.3777  E-5  mm.  correspondiAndole  una  frecuencia 
predominante  de  0.719  hertz,  la  amplitud  minima  obtenida 
del  analisis  espectral  fue  de  0.  421^9  E-6  mm.  sicndo  r.u 

frecuencia  de  1.344  hertz. 

Su  frecuencia  de  corte  cor respondiente  es  de  1.039  hertz. 
Se  puede  observar  en  la  figura  12.8  Csismo  62),  qus  la 
amplitud  mAxima  del  espectro  fue  de  0.4861  E-5  mm.  y  la 
frecuencia  predominante  que  le  corresponde  de  0.813 
hertz,  la  amplitud  minima  fue  de  0.4460  E-6  mm.  y  su 

frecuencia  correspondiente  de  1.219  hertz. 

La  frecuencia  de  corte  obtenida  del  espectro  es  de  1.01 
hertz . 

El  espectro  de  la  figura  12.9  Csismo  70)  tambicn  nos 
muestra  una  amplitud  mAxima  del  espectro  de  0.5783  E-5  mm. 
y  su  frecuencia  es  de  0.844  hertz,  la  amplitud  minima  cs 
de  0,4573  E-6  mm.,  siendo  su  frecuencia  de  1.344  hertz. 

La  frecuencia  de  corte  correspondiente  es  de  1.024  hertz. 
La  figura  12.10  Csismo  64)  tambiAn  nos  muestra  unn 
amplitud  mAxima  y  minima  del  espectro  de  0.3713  E-5  y 

0,5965  E-6  mm.  y  sus  correspond ientes  frecuencias 

predominantes  de  0.306  y  1.313  hertz  respect ivamcntc. 

La  frecuencia  de  corte  obtenida  es  de  1.054  hertz. 

El  espectro  12.11  Csismo  1 16)  tiene  una  amplitud  mAxima  do 
0.6475  E-5  mm.  y  su  frecuencia  predominante  de  0.750 
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cd  m  io 


hertz,  siendo  la  minima  amplitud  de  0.2723  E-6  mm.  de 

frecuencia  1.25  hertz. 

La  frecuencia  de  corte  os  de  0.39  hertz. 

-  La  figura  12.12  (si smo  118)  tiene  la  amplitud  maxima  de 
0.9419  E-5  mm.  y  su  frecuencia  predominante 

correspondiente  de  0.906  hertz,  siendo  la  minima  amplitud 
de  0.9349  E-6  mm.  y  su  frecuencia  1.188  hertz. 

Su  frecuencia  d c  corte  es  de  0.969  hertz. 

Adem£s  los  espectros  de  la  segunda  parte  muestran  una  fuerte 
caida  de  las  amplitudes  que  es  similar  a  los  espectros  de  la 
primera  parte  analizada. 

TABLA  N2  9 


AMPLITUD  FRECUENCIA  FRECUENCIA  AMPLITUD  FRECUENCIA  VEL. 
MAXIMA  <AMPL.  MAX.)  DE  CORTE  MINIMA  CAMPL.  MIN.)  SPEC. 


SI  SMO 

E-05  mm. 

hertz 

hertz 

E-06  mm. 

hertz 

Km/s 

7 

0.4214 

0.688 

1 . 039 

0.1117 

1.  156 

3.  15 

12 

0.2097 

0.750 

0.773 

0.3171 

1.219 

3.  18 

16 

0.9714 

0.688 

0.906 

0.4335 

0.969 

3.  18 

22 

0.8545 

0.898 

0.930 

0. 1457 

1.365 

3.06 

62 

0.4861 

0.812 

1.012 

0. 4-160 

1.219 

3.  16 

64 

0.9719 

0.906 

1.054 

0.5965 

1.313 

3.  17 

70 

0.5783 

0.844 

1.024 

0.4573 

1 . 344 

3.  10 

72 

0.9777 

0.719 

1.039 

0.4218 

1.344 

3.  IS 

78 

0.5275 

0.719 

0.990 

0.  1643 

0.906 

3.  12 

93 

0. 1844 

0.742 

0.915 

0. 1852 

1.055 

O  •  X  4tm 

116 

0.6474  ' 

0.750 

0.890 

0.2723 

1 . 250 

3.03 

118 

0.9419 

0.906 

0.969 

0.9349 

1.  188 

3.  19 

En 

general 

podemos 

decir  que  las 

car act er 

i st  icas 

de  los 

espectros  obtenidos  nos  muestran  que: 

-  Las  met»imas  amplitudes  espectrales  se  encuontran  dentro 
del  rango  de  0.1844  y  0.9777  E-5  mm. 

-  Quo  las  frecuencias  pr edominantes  (corr espond ientes  a  las 
amplitudes  mAximas)  se  encuentran  dentro  del  rango  de 
0.688  y  0.906  hertz. 
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Las  amplitudes  espectrales  minimas  est&n  contenidas  cntrc 
los  0.1117  y  0. 33-49  E— 6  mm. 

Las  frecuencias  predominantes  (do  las  amplitudes  minimas  > 
son  de  0.906  a  1.344  hertz. 

Las  frecuencias  de  corte  encontradas  estAn  en  el  range  de 
0.773  a  1,059  hertz.  ^ 

La  parte  plana  de  los  espectros  estA  mostrada  por  la  linea 
trazada  sobre  las  amplitudes  m&s  predominantes  del 
espectro  fig.  12.1  al  12.12,  y  los  valores  obtenidos  del 
grAfico  estAn  dentro  del  rango  de  0.3  E-5  a  0.9  E-5  mm. 

Se  puede  observar  en  cada  figura,  que  la  parte  plana  del 
espectro  tiene  un  contenido  menor  de  amplitudes,  cn 
relacidn  con  la  parte  obi icua  que  contiene  las  amplitudes 
que  decrecen;  ello  nos  indicar  la  que  la  energia  que  se 
mantiene  en  la  parte  plana,  es  menor  en  relaciAn  con  la 
energia  que  se  pierde,  lo  cual  estA  most r ado  en  la  parte 
obi icua  del  espectro. 

En  la  intersecc iiin  de  la  parte  plana  y  la  parte  obi  icua  de 
los  espectros,  tenemos  la  frecuencia  de  corte,  que  es  la 
que  determina  el  punto  a  partir  del  cual  decrees  la 
energia. 

Sc  ha  observado  que  las  frecuencias  predominantes  que  se 
obtuvieron  de  los  espectros,  tanto  para  sismos  de 
recorrido  continental  y  el  sismo  que  se  encuentra  proximo 
a  la  linea  de  costa  (recorrido  oceAnico),  no  muestran 
diferencia  alguna,  lo  que  nos  indie aria  que  el  sismo  quo 
sc?  encuentra  cerca  a  la  linea  do  costa  tiene  su  recorrido 
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a  travel 
si  mi  1  arcs 
onuuontren 
Cost  era  y 


Fig. 


<|r>  cr.tructuran  con  prop  i  odados  olAsticas 
a  los  sismos  de  recorr ido  continental,  quo  sc 
on  las  provincias  geol  <‘g  i  c  as  do  la  LLanura 
Cor  d ill cr a  Oc c i don  t  al  . 
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11.-  Loc  al i 2 ac i 6n  do  espec t r os  anal i z  ados . 
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Cap i tul o  VI 


Posible  Mecanismo  de  Tr  anninisiiin  de  1  as  on  das  Rq 

V 1 . 1 . -  F und ament  os  tear Icos 

Para  tener  una  vision  cor  recta  para  L_,ca  parte;  del 
trabajo,  os  nscesario  introduc  ir  ^algunos  f undamen tos 
teir  icos,  qu.e  a  la  vez  nos  servirAn  para  captar  lex  idea 
del  posible  mecanismo  de  tr ansmisiAn  de  la  Rgf  cuyo 
ent r ampamiento  se  prodaciria  en  una  capa  quia  dentro  de 
la  corteea  terrestre;  para  ello  tendremos  a  las  : 

V I . 1 . 1 . -  Ondas  Guiadas 

Se  llaman  as!  las  per turbac  Lories  que  se  transmiten 
trav^s  de  un  cuerpo  elAstico  y  se  prcpaga.i  en  medico 
formados  de  estratos  hor isontal es,  en  los  quo  sc 
notan  cambios  de  propiedadss  en  sent i do  vertical,  ma¬ 
ne  en  el  hor iron led.  Las  ondas  super f i c i al es  son  ondas 
guiadas  en  la  cortesa  terrestre  y  a  este  tipo  de  ondas 
pertenecen  las  ondas  Love  y  Rayleigh. 

VI. 1.2.-  Ondas  Canal es 

Ocurre  que  en  ciertos  ni voles  la  eneryia  de  las  ondas 
se  c one ent r a  transmit iAndoso  mejor  que  en  n i vel os 
adyacentes,  deb i do  a  las  pr opiedades  c?l  Ast icas  del 
medio;  las  prov.  iinidades  de  este  ni vel  critics 

const  it uyen  un  canal  y  las  ondas  guiadas  que  so 
transmiten  por  Al  son  las  ondas  canal  es  CPayo,  1  Of  0  *  . 
Una  capa  de  ba.ja  velocidad,  que  cst  A  const  i  tui.J  \  p ...  r  un 
medio  de  distintas  >_•  -ns  t  antes  cl  Ast  icas  pued;. 
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comport arse  como  un  canal  por  el  cual  se  transmiten 
ciertas  ondas  canal  es  y  den  t  r  o  cJe  este  tipo  do  oridati 
guiadas  podemos  mencionar  las  siguientes: 

La  F'g  y  Sg  son  ondas  que  r.e  propagan  cn  la  capa 
gr an i  t  i  c  a . 

Las  Px  y  Oh  r,on  ondas  quo  r.e  pr*:<pagan  on  la  capa 
basAl t ica. 

Las  Pn  y  Sn  son  P  y  S  que  se  propagan  por  un  canal  cn 
la  parte  superior  del  man  to. 

Las  ondas  T  son  ondas  guiadas  por  uri  canal  dc  ba.ja 
velocidad  que  es  el  oc6ano. 

Las  Li  son  ondas  que  se  propagan  on  la  capa  dc  baja 
vel oc i dad  del  raanto  super i or . 

Las  Lg  son  ondas  super f i c i al es  dc  cor  to  per iodo 
que  se  propagan  por  la  capa  gran  it ica. 

VI. 1.3.-  Ref 1  ex i ones  Mtil t ipl es 

Las  ondas  sismicas  originadas  en  el  foco  son  on  parte 
ref lejadas  a  la  super ficie  desde  el  contacto  do  dos 
capas  de  velocidades  diferentes.  Segijn  el  principio  de 
Huygens,  cada  punto  del  contacto  alcansado  por  la  onda 
puede  considerarse  como  un  nuevo  foco  dc  energia  que  tla 
origen  a  nuevas  ondas  elAsticas,  que  se  propagan  segi'm 
las  1  eyes  de  reflexion  y  refract  i<in,  de  la  misma  manera 
que  la  onda  primitiva. 

Como  siempre  existe  el  contr aste  de  velocidades 
entire  la  capa  mAs  consol  idada  (subyacente  t  y  la  capa 
mAs  superficial  ' ::upr  ayacentt?)  ,  la  onda  originada  en  el 
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contacto  puede  ser  otra  ve:  reflejada  desde  la 

super ficie  o  desde  la  capa  mAs  superficial  hacia  abajo. 
Dicho  proccso  puede  repet irso  varias  voces  on 

condic  iones  favorables  y  a  estas  reflexionoa  sc  lot. 
llama  ref lex iones  multiples.  En  la  figura  (1C)  tenonos 
quo  c  os  cl  receptor  do  las  ui.das  yz  cc  el  foco  del 
sismo;  supongamos  quo  la  vclocidad  es  const ante  <  B*  ) 
on  la  capa  guia,  la  onda  reflejada  sigue  el  camino 
sin  embargo  es  posible  que  se  produzca  la  reflexion 
multiple  quo  seguirA  el  camino  adefc  que  sc  registrar  A 
en  c  con  el  tiempo  doble  que  la  abc. 

Desde  luego  que  la  reflexiAn  multiple  adefc  no 
la  (mica  posible;  si  el  reflector  es  bueno  y  la  energia 
suficiente  se  producirA  la  reflexiAn  mdl t ipl cahkl emnpc 
y  aCin  otras;  de  esta  manera  las  reflexiones  multiples 
aparocerAn  per i Adicamente  y  podrAn  ser  i  dent  i  f  i  c  adas 


Fig.  13.--  Re  f  1 x  iones  multiples 
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VI.  1.4 


Inter  ferenc i a  Construct Iva 


Supongamos  una  capa  du  espesor  It  y  por  dobajo  un  medio 
semi— in f i n i to  siendo  las  velocidades  un  las  capas  < 
Cfigura  14.>  y  consider ando  que  los  rayos  quo  inciden 
ert  la  super  ficie  de  separ ac  i An  tienen  un  Angulo  de 
cmergenc  ia  menor  que  el  Angulo  de  enter  gene  i  a  cr  it  ico 
Ce  <  e’>,  consider omos  un  punto  en  la  capa,  tal  que  en 
un  memento  dado  su  movimiento  serA  la  suma  de  todos  los 
movimientos  de  las  ondas  que  pasen  por  £1 .  Si  varias 
01  id  as  estAn  en  fase,  la  inter  Ferenc  ia  amplia  la  energia 
del  movimiento  y  se  cl.cc  que  es  construct  iva,  en  cambio 
si  estAn  en  oposiciAn  de  fase  disminuye  la  ampl itud  y 
se  dice  que  la  inter ferenc ia  es  destr uct iva. 


Su|«mI  !«•  I  »•  I  litre  I  ttl  »*r  IVri'ite  la  Const  rurl  I  vii 


Fig.  14.-  Inter  ferenc  ia  construct  iva. 

V I . 2 . -  Fosibl e  Mecan i smo  de  Tr ansmisi An 

Son  pocos  los  trabajos  real i 2 ados  sob re  mecanicmo  de 
transmisi An  de  este  tipo  de  ondas,  sin  embargo  Path 
(1954)  nos  indica  quo  las  ondas  Rg  existen  para  c ami nos 
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enter amente  cont inental os,  quo  la  distr  ibuc 1  An  de  --u., 
velocidades  es  independiente  do  la  distancia  opicentrnl 
y  que  se  propagan  on  una  capa  de  baja  vclocidad,  quo 
s e  encuentra  eri  la  corteza  terrestre,pero  su  estudio  io 
real iza  on  base  a  o ndas  quo  existen  para  profundidades 
mayor  os  a  la  normal;  ougicrc  quo  existc  on  la  cor  to.:  a 
terrestre  una  capa  gu i a ,  que  es  la  responsable  de  1  a 
transmisiAn  de  este  tipo  de  ondas. 

« 

Gutenberg  (1351)  infiere  la  existencia  de  una  capa 
de  ba.ja  velocidad  a  una  profundidad  de  15  km.,  debido  a 
discrepancies  entre  las  velocidades  de  las  ondas 
longitudinales  obtenidas  de  terremotos  corcanos  y 
explosiones.  Tambi6n  entre  el  anAlisis  de  las  ondas 
longitudinales  y  transversal  os  r eg i sir  ados  de 
terremotos,  cn  los  cuales  ocurre  un  deer ec  imi onto  do 
las  const ant os  elastic  as  con  la  profundidad  dobido  al 
incremento  de  la  temper atura.  A  estas  ondas  do 
velocidad  constanto  que  so  transmiten  por  la  capa  de 
ba.ja  velocidad  las  denomina  ondas  canales.  Para 
nuestro  model  o  de  mocanismo  do  pr  opagac  i  An  do  las  on  dev.. 
Rg  se,  deben  explicar  las  car  actor  1st  icas  y  prop  icidados 
que  se  ban  encontr ado  en  capitulos  an  tor iores  y  aci 
esbozar  un  posible  modclo  do  pr opagac i An  do  las  ondas 
en  estudio;  por  tan to  hay  quo  explicar  : 

1 .  -Cub  la  onda  Rg  se  defina  mo.jor  on  c  ami  nos 

c ont inentales,  mas  no  on  c ami nos  oceanic os. 

2.  -Due  ol  contenido  do  amplitudes  del  tren  do  ondas  do 
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la  Rg  tonga  una  duraci<5n  on  tiempo  do  213  s . 

3.  ~E1  hccho  quo  scan  on  das  elipticas  y  rctr  <5qr  acJas. 

4. -Explicar  el  hecho  que  sean  ondac  do  ba.ja  volocidad. 

5.  -Que  las  frccucncias  varlcn  y  quo  sus  amplitudes  m 

vayan  atenuando  a  partir  do  la  tercera  zona. 

6.  -El  hecho  quo  on  los  espectros  1  £\  parte  plana  sea 

monor  que  la  parte  oblicua. 

-  Para  la  region  do  estudio  ha  sido  posible  registrar 

* 

sismos  con  fases  Rg  de  recorrido  oceAnico  y 
continental;  por  ello  so  analizaruri  por  sepav  ado.  Do 
los  sismos  de  recorrido  oceAnico  se  ha  podido  observar 
para  la  primer a  zona  que: 

La  fast?  Rg  se  observa  para  epicentres  en  el  Pacifico; 
para  distancias  de  0,2°  a  1,62  exist iondo  gran 
atenuaciAn  con  excepci.fm  del  sismo  05,  que  prcscnta  una 
fase  de  Rg  clara  (de  0,552  de  distancia). 

Do  la  linea  de  costa  hasta  una  distancia  do  0,22,  el 
desarrol lo  de  Rg  se  presents  con  poca  atenuaciin. 
AdemAs  es  necesario  consider ar  quo  para  la  sogunda  y 
tercera  zona  eri  todos  los  sismos  ubicados  en  el  ocPano 
se  observ-5  alta  atenuacitfn  do  Rg,  espec  i al mente 
en  aquellos  sismos  cuyos  focos  se  encuentran  en  la 
tercera  zona. 

El  model o  expuesto  por  Fisher  y  Raitt  (1962),  menciona 
quo  el  espesor  de  la  corteza  vnria  a  medida  quo  so 
alcja  do  la  linea  de  costa  al  ocPano  (17  a  11  Ini.), 
s iondo  variable  en  toda  la  linen  do  costa  on  estudio; 
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si  consider amos  ello  podromcc  espl  icar  l.v: 
car  acter  i  st  icas  raencionadas  art  ter  ior  monte,  de  l.i 
siguiente  manera: 

Sabemos  quo  los  mejores  registros  se  obtuvicron  cerca 
de  la  linea  de  costa,  mas  no  los  que  se  cncuentran  a 
mayor  distancia  de  0,29,  se  podria  exp 1 i c ar  por  l as 
diferencias  de  espesores  existentes  entre  la  linea  do 

costa  al  oedano  de  la  region  dc  estudio.  Giendo  los 

4 

sismos  mds  prrtximos  a  la  costa  los  que  mejor  so 
canal  i^en  on  lex  capa  guia  dc  la  cortcna  y  tengan  mc'nor 
pdrdida  de  energia.  Si  consideramos  esto  exp  1 icar i amos 
el  hecho  de  que  las  ondas  a  mayor  distancia  do  la  linea 
scan  sus  llegadas  un  tanto  atenuadas  debido  a  que  su 
canal isac itfn  habria  sucedido  cn  la  corteza  oceAnica, 
prAxima  a  la  linea  de  costa  y  en  su  recorrido  a  travds 
de  el  la  perderia  energia  que  seria  absorb ida  por  el 
manto  y  el  oedano  (fig.  15. a  ,  15. b). 

Para  los  de  recorrido  continental  de  la  primera  y 
segunda  zona  se  obtuvier >f«n  buenos  registros  de  la  faso 
Rg,  aun  cuando  existen  algunos  sismos  atenuados  (sismos 
56  y  65)  ubicados  en  la  Cordillera  y  Al t i piano.  Csto 
puede  debersc  al  paso  de  las  ondas  Rg  por  debajo  de  la 
Cordillera  de  los  Andes,  donde  la  corteza  nos  muestra 
una  estructura  irregular.  Las  ondas  de  1  u  t  o  v  c  c?  v  a  r  o  1 1  ex 
present an  mayor  recorrido  continental  y  alta  atcnuaciin 
pociblemente  debido  al  paso  do  las  ondas  Rg  por  la 
corteza  irregular.  Si  consider amos  lo  oxpuesto,  la  Rg 
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ser ia  una  face  cuya  capa  guia  que  canal ina  lac  ondac 
estaria  c:n  la  cor  tec  a. 

-  Que  el  contcnido  de  amplitudes  del  tren  do  ondas  Rg 
tenga  una  durariAn  de  20  s  indica  que  dcbe  exist ir 
una  capa  guia  de  baja  velocidad  que  sea  la  cauoante  de 
su  canal  izac  i  >5n,  ocas  ion  an  do  que  haya  poca  p^rdida  de 
energia  durante  su  transmisitin  del  foco  a  la  estaciin 
de  registro. 

-  El  movimiento  de  particula  de*  la  Rg  ticne  un 
movimi onto  eliptico  retrograde,  el  cual  nos  indica  que 
son  ondac  de  tipo  Rayleigh,  por  lo  cual  son  ondac 
super  fic  iales  y  que  los  causantes  de  su  formaciOri 
seriari  las  ondac  de  cuter po;  s iempre  que  cstas  ondac 
corpOreas  estGn  confinadas  al  limite  super  fide  libre- 
corteza  y  las  ondac  so  polar  icon  de  modo  que  las 
particulas  del  medio  se  rouevan  en  el  piano  vertical 
paralelo  a  la  dirccciAn  del  movimiento  de  la  onda. 

-  La  baja  velocidad  de  Rg,  indica  que  son  ondac  cuya 
transmisi  6n  se  hace  a  travAs  de  una  apa  superficial 
donde  puedan  transmit ir so  con  poca  pdrdidn  de  energia. 
-Si  consideramos  a  la  cortez a  como  medio  de  grandee 
irr egul ar idades  expuesta  a  las  zonas  do  baja  velocidad 
de  la  capa  granitica  (Qcola,  1972),  a  1«  complejidad  do 
las  raiceo  de  los  Andes,  al  inagmat  ismo  producto  del 
vulcanismo  existontc.-,  todo  olio  nos  podria  dar  una 

expl icasi An  de  la  var iaciAn  de  las  f rccucnc ins  y  la 
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at cnuac i in  de  las  ondas  en  la  tercera  zona,  ya  quc  !a 
onda  Rg  estA  atravezarido  un  medio  complejo  de  grander 
i nhomogene i d  ades . 

-  Este  hecho  nos  muestra,  como  la  energia  que  sc  pierde 
es  mayor  quo  la  energia  quo  se  mant iene (par te  plana), 
csta  per  did. a  do  energia  sc  lc  atribuya  al  recorr  id... 
por  debajo  do  la  Cordillera  de  los  Andes. 

Si  consider amos  estas  car ac ter  1st icas,  dir emos  quc  la 
onda  Rg  ha  sido  produeto  de  las  ondas  de  cucrpoj  quo  al 
incidir  la  onda  longitudinal  en  una  discont inuidad 
entre  dos  medios  (super  ficic  libre  --  cor  to.: a. > 
producirian  ondas  P  y  SV  ,  quo  se  canal izari an  en  la 
cortez  a  lie  gando  a  la  est  ac  i  >f.n  de  registro  por 

constantes  reflex iones  y  ref race ioncs;  para  quc  se 
mantenga  la  energia  observada  en  la  parte  plana  de  1 
espectros  sc  tendria  quc  produeirsc  la  inter  ferciv:  in 
construct iva  en  la  capa  guia  que  es  el  medio  por  el  quc 
se  transmiten  las  ondas  Rg. 

De  acuerdo  a  lo  anterior,  nos  faltaria  determinar  a  que 
profundidad  estaria  propagandosc  las  ondas  Rg;  para 
ello  asumiremos  el  model o  de  corteza  de  James  (1971) 
como  tambiC'ii  de  Ucola  y  Meyer  (1972),  para  la  pm  te  del 
altiplano,  deba  jo  de  Per i_V--Bol  i via. 

En  esta  tabla  nos  muestra  las  velocidadcs  de  los 
primeros  arribos;  ahora  bien  consider ando  la  rclaci4n  do 
velocidadcs  de  las  ondas  principales  como  cl  radio  cir 
poisson  de  1  /A  se  put.  de  ob  tenor  la  vcl  oc  i  dad  do  la  onda 
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para  ambos  per  files  (a  y  b>  con  la  siguionte 


rel  ac  i  -in: 


Vs  ~  Vp  /  3 


) 


TABLA  #  9 


(model u  do 

c  or  t  ez  a 

para  ol 

A1  tipi  a no 

sur  dot  Por  Ou 

DEPTH 

Va 

Vb 

Vs  <  a ) 

Vs  ( b  j 

km. 

k  m  /  s 

km/s 

km/s 

km/s 

0 

4.50 

4.50 

2.59 

2.59 

6.4 

4.50 

4.50 

2.59. 

2.59 

6.5 

6.04 

6.01 

3.49 

3.47 

9.0 

6.04 

6.04 

3.49 

3.49 

9.  1 

5.00 

5.00 

2.89 

2.89 

11.0 

5.00 

5.00 

2.09 

2  ■  S3 

12.0 

6.  10 

6.08 

'"i  rr 

vJ  • 

3.51 

30.5 

6.  10 

6.  12 

3.52 

r? 

s~t  m 

31. 5 

6.  75 

6.75 

3.  90 

3.90 

36.5 

6.75 

6.70 

3.90 

3.91 

36.6 

6.  15 

6.  15 
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6.  15 

6.  15 

3.55 

3-  55 

45.5 

6.90 

6.90 

3.98 

3.98 

72.0 

6.90 

6.90 

3.  98 

3.  98 

73.0 

0.00 

8.00 

4.62 

4 . 62 

90.  a 

8.  00 

8.00 

4.62 

4 . 62 

Si  consider amos  las  velocidades  de  las  ondas  Rg  ;  el  las 
varian  entre  los  2.05  —  3.40  km/s  y  las  relac ionamos  con 
las  veloci  (Jades  de  los  model  os  de  James,  Ocola  y  Meyer 
se  puede  observar  que  se  encucntran  dentro  del  rango  de 
velocidades  de  2. 59  -  3.49  obtenidas  por  Ocola  y  Meyer 
como  tambien  dentro  del  rango  de  velocidades  obtenidas 
por  James  2.80  -  3.4C  km/s.  Elio  hace  sup oner  quo  la  Ry 
c5  probable  que  se  cncuentre  propagandose  por  una  capa 
quia  de  2.6  km  de  ospcsor  la  cual  so  encontraria  por 
cncima  de  la  capa  de  tn.ja  vrlocidad  pair  a  el  model  o  do 
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Ocola  (fiqma  16. a),  y  tin  osposor  rfc  G.U  l:m  para  t?.l 
model o  tie  Jamcc  (tablaG  It  1  ,.,1  7),  amt)  os  dentro  de  la 
cap a  granitica  (fig.  16. b). 
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FOC  0 


l.ineu  tie 


Costa 


HO  C  0 
t- . 


p  j  q .  J  „  3  >  F  oc o  & I  *.?  ttc.l tie*  1  3  l  do*  ^ o  • 

b)  Face  cerca  d«j  Jo  l  inos  do  Co-,to . 
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6.4  km . 


Capa  Sedimentaria 


Capa  Guia  (posible  propagacion  de  Rg)  j 

2.6 

2.0 

Capa  de  baia  velocidad  j 

8.0 

Capa  de  baja  velocidad 

5.0 

(a) 


Capa  Sedimentaria  * 

5 . 0  km . 

Capa  Guxa  (posible  propagacion  de  Rg) 

6 . 0  km . 

6.0  km. 

( 

8 . 0  km . 

8.0  km. 

(b) 

Fie.  16.-  Modelos  tornados  para  tin  posible  mecnmismo  dc  transmlsion 
de  Rg. 

(a)  Modelo  tornado  de  Ocola  1972  (1)  Capa  granitica  (2) 
Capa  Gabroide 

(b)  Modelo  tomadp  de  James,  1971  (l)  Capa  granitica  (2) 
Capa  Cabroide 


Cone  1  us i ones  y  Recomendac iones 

Salvo  nxcopc  ionrs  rl  comirnzo  or,  eracryonto. 

-  Los  periodos  van  do  0.0  a  1.5  segundos. 

Las  veloc idades  varian  entre  £.8  a  3.4  km/s. 

-  Las  mayores  amplitudes  se  observa  en  la  zona  sur. 

~  Las  pro fund idades  ostAn  outre  super f ic iales  e  intermedins 

-  Las  amplitudes  de  los  sismos  con  epicentres  en  el  Delano  para 
la  zona  sur  es  variada,  siendo  de  mayor  predominant:  ia  para 
sismos  que  se  encuentran  cere  a  a  la  Costa. 

-  Las  amplitudes  de  sismos  con  epicentres  en  el  Oc^ano,  para  la 
zona  central  y  norte  muestran  alta  atenuac iAn. 

-  No  se  adviorten  diferencias  en  las  diversas  provincias 

geolAgicas  dentro  de  la  zona  sur.  * 

-  En  la  zona  central  hay  diferencias  mar  cad  as  sogCin  las 
provincias  geolAgicas;  las  mayores  amplitudes  en  sismos  del 
Oubandino,  las  mlnimas  amplitudes  en  sismos  do  la  Cordillera 
Occidental  y  linea  de  Costa. 

-  La  zona  norte  es  la  presents  mayor  comple.jidad  en  relaciAn  a  la 
amplitud  de  Rg;  aqul  se  observe  amplitudes  menores  a  las  otras 
d  o  s  z  o  n  a  s . 

-  El  anAlisis  del  movimiento  de  particula  presents 
irregular idades,  pero  prevalece  el  movimiento  vertical 
longitudinal  y  se  dibuja  el  movimiento  eliptico  retrograde 
caracter  1st  ico  de  las  ondas  de  Rayleigh. 

-  La  polar izaciAn  aparente  est A  comprendida  entre  3.52  y  422 
respecto  del  aziim.it  Estacion—  Epicentre. 

-  L:v.,  mAximas  amplitudes  espectrales  estAn  comprendida::,  entre 
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-  Se  interprets  que  las  ondas  Rg  se  propagan  por  la  Cortez  a,  [3  or 
una  capa  guia  que  es  probable  se  encuentre  en  la  capa 
granitica  y  est£  propagandose  por  una  capa  de  espesor  de  3  km 
de  espesor. 

-  El  subsuelo  marine  atenua  notablemente  las  ondas  Rg,  como 
tamb i6n  la  complejidad  del  subsuelo  Cordil lerano  que  present n 
atenuaci'in  no  tan  fuerte  como  el  subsuslc*  marine. 

RECOMENDAC I ONES 


-  Es  recomendable  extender  el  presente  estudio  a  toda  Sudambrica. 

-  Convendr ia  hacer  estudios  sobre  la  coda  del  sismo  para  hacer 
coinparac  iones  con  los  t rones  de  ondas  de  la  Fig. 

-  Ademas,  estudios  de  dispersion  de  ondas  Rg  para  oncontrar  un 
model o  de  corteza  mejorado;  por  el  cual  se  propaguen  este  tipo 
de  ondas. 
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ANEXO  N°  1 


G  » 

10  REM  ****  ****************  ******************  **  ******  **  *  *******  k  **»*  ***  * 
20  REM  *********  PROGRAMA  PARA  GRAF  ICAR  MOVIMIENTO  DE  ******  >•  ************ 

20  REM  *******  PARI  ICLJL A  DR  LAG  ONI) AG  Lg  *****  ************ *  * 

40  REM  *******  22 -XI  —GO  Rodolfo  Ayala  *****  *:i  x  ********* ** 

GW  REM  *** *  ***************  **********  *******  *********  *****  ****  * *** ** *  *  *  * 
GW  REN 

70  DIM  AC  1200  ) , BC 1200) , NX  19001 , NY  C900) 

30  D I M  C  C 1 200 > , PH < SOW > , PV  C 900  > , MX  t  C 90 > , NY l C 90 ) 

90  CTO 

100  PRINT  "  M 

110  INPUT  "NOMBBE  DE  ARCHIVO  :  "  ,  X  T> 

120  PRINT  "  " 

130  INPUT  "  NUMERO  TOTAL  DE  PIJNTOS  :  " ,  N 

140  INPUT  "  NUMERO  DE  GISMOS  :  " , D 

100  PRINT  "  " 

1G0  INPUT  "  A2IMUT  CEpi-Est)  :  " ,  A7.T 

170  CLS 

100  LOCATE  11,25 

190  PRINT  " ! ESPERAR !  ADQUIR IENDQ  DATOS  " 

200  OPEN  "I"  ,411,  X* 

210  INPUT  0 1 , A$ 

220  FOR  I-  1  TO  N 
230  INPUT  111,  AC  I) 

240  NEXT  I 
250  N-N/3 
2 GO  S--AC1) 

270  FOR  I “2  TO  N 
'i'80  S=A  C I  .>  +5 
290  NEXT  I 
300  0=5/ N 
310  FOR  1=1  TO  N 
320  aci:>=aci:>-5 

330  NEXT  I  a 

340  J=0 

350  FOR  I-N+l  TO  2*N 
3G0  J=J  H 
370  BCJi-ACI) 

300  NEXT  I 
390  S=B  C 1 ) 

400  FOR  1=2  TO  N 
4 1 0  0=D C I > +3 
420  NEXT  I 
430  5=5 /N 
440  FOR  1=1  TO  N 
450  B  C I  .)  =D  C I  >  -S 
4G0  NEXT  I 
470  |<=0 

430  FOR  I=2*N+1  TO  2*N 
490  F=K+1 
500  C  CK>  =  -A  <  I  > 

510  'PRINT  CCL i , CC 1 1 
520  NEXT  I 
530  3  -r  C 1  ;> 

5-10  FDR  T  '•  TO  1 1 
550  5  -C  '  t  '  :  0 
5C0  NEXT  I 
5/0  5  5/N 
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580  FOR  1=1  TO  N 
580  C ( I >  =C C I > ~S 
600  NEXT  I 

610  IF  AZT:>  .100  THEN  AZU-AZTH80 
620  IF  AZT>  100  THEN  AZU=AZT-1O0 
600  AZI  =  <AZU*3. 1 1 1 529 >  X 1 00 
640  FOR  1=  1  TO  N 

650  PH  <  I  >  =B  (  I  >  *GQS  <  AZ I )  +C  C  I  >  x  3 1 N  C  A 1 1  ) 

660  PV  ( I  )  =  D  <  I  :>  *  S INCAZI )  --C  C 1 >  * COS  (AZ  I  > 

670  NEXT  I 
600  CLS 

690  INPUT  "CONTINUAR  (0)  ;  IMPRIMIR  Cl)  ;  TERM 1 MAR  (2)  :  X 

700  IF  X=0  THEN  GOTO  790 

710  IF  X  - 1  THEN  LOCATE  21,1:  PR  I  NT  "  OPR  I  MIR  :  1MPR  PANT'  " 

720  'WAIT 
730  LOCATE  21,1 
740  PRINT  " 

750  LOCATE  22, 1 
760  PRINT  " 

770  GOTO  2090 

780  IF  X=2  THEN  GOTO  2090 
790  PRINT  "  " 

800  INPUT  "PUNTO  INICIAL  Y  FINAL  :  ",Y,N 

oi0  amax--a<y:> 

820  AM I N= AMAX 
830  FOR  I=Y  TO  N 
840  T = A ( I > 

850  IF (T  AMIN) <0  THEN  GOTO  830 
860  IF  (T-AMIN)  >=0  THEN  GOTO  090 
870  IF  (T-AMIN) >=0  THEN  GOTO  900 
880  AMIN--T 

090  IF  CT-AMAX ) <  =0  THEN  GOTO  920 
900  IF  (T  AM AX)  0  THEN  GOTO  910 

910  AMAX--T  « 

920  NEXT  I 
930  AM AX  1=0 (Y ) 

940  AMI  HI “  Ah AX  1 
950  FOR  I =Y  TO  N 
960  T 1=D ( I ) 

970  IF  ( T 1  --AM I N 1  )  < 0  THEN  GOTO  1000 
980  IF ( T 1  AMI N 1 ) >=0  THEN  GOTO  1010 
990  I F ( T 1 -AM  I N l > >=0  THEN  GOTO  1 020 
1000  AM I N 1 =T 1 

1010  IF  ( T 1  --AMAX 1 )  <  =0  THEN  GOTO  1040 
1020  IF  CT1-AMAX1) >0  THEN  GOTO  1030 
1030  AMAX 1 =T 1 
1040  NEXT  I 
1050  AMAX2~:C  <  Y  ) 

1060  AM I  M2 "AMAX 2 
1070  FOR  I-Y  TO  N 
1000  T7  At  f  > 

1090  IF(T2  AMIN2) <0  THEN  GOTO  1120 
11.00  IF  (12-  AMIN2  >  0  THEN  GOTO  1130 
1110  IF  (  T2  AM  I N2  >  >r-0  THEN  GOTO  114  0 
1.120  AM  I  M2- *T  2 

1130  IF  (T 2  -AMAX 2)  '■  =0  THEM  GOTO  1160 
1140  IF  (T2-  AMAX 2)  >0  THEN  GOTO  1160 
1150  AMAX  -T 
1.160  NEXT  I 
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1170  AMAX3=PH(Y) 

1100  AMIM3=AMAX3 
1190  FOR  I =Y  TO  N 
1200  T3=PH C I ) 

1210  IF (T3-AMIN3) <0  THEN  GOTO  1240 
1220  IF  (T3-AMIN3)  >=0  .THEN  GOTO  125)0 
1230  IF  (T3--AMIM3)  >--0  THEN  GOTO  1260 
1240  AM I N3”T3 

1250  IF  < T3—AMAX  3) < *0  THEN  GOTO  1280 
1260  IF  (T3-AMAX3)  >0  THEN  GOTO  1270 
1270  AMAX3--T3 
1200  NEXT  I 
1290  AMAX4=PV(Y) 

1300  AMTN4--AHAX4 
131.0  FOR  I  —  Y  TO  N 
1320  T4=PVCI) 

1330  I F ( T 4--AMIN4) <0  THEN  GOTO  1360 
1340  IF  (T4-AMIN4) >=0  THEN  GOTO  1370 
1350  IF  (T4-AMIN4) >~0  THEN  GOTO  1330 
1360  AMIN4-T4 

1370  IF  CT4--AMAX4)  <~0  THEN  GOTO  1400 
1380  IF  (T4-AMAX4 ) >0  THEM  GOTO  1390 
1390  AMAX4=T4 
1400  NEXT  I 
1410  PRINT  "  " 

1420  PRINT  SPC(1)"Z"  SPC( 13)  "N"  5PCC13)  *'E"  SPC(IG)  "PH"  SPG  Cl 
1430  ’  FOR  I-Y  TO  N 
1440  PRINT  "  " 

1450  PRINT  AMAX, AMAX1, AMAX2, AMAX3, AMAX4 

1460  PRINT  AMIN, AM INI, AMI N2, AMI N3, AMI N4 

1470  ’ OPTION  DASE  1 

1430  KEY  OFF 

1490  SCREEN  2 

1500  PRINT  "  " 

1510  INPUT  "VALOR  MAXIMO  EBCAl.A  (POSIT  I  VO  0  NEGAHvb)  :  ",  MAX 
1520  MAXZ~  AB3CMAX) 

1530  PRINT  "  " 

1540  INPUT  "VALOR  MARC A  DEL  IN I CIO  (1-5)  :  ",Q 

1550  CLS  0 

1560  LOCATE  1,15 

1570  PRINT  "SISMO  No.  :  "D 

1500  LOCATE  1,45 

1590  PRINT  "A? (Est-Epi )  :  "AZU 

1600  REM  **********  PRO  YECC I  ON  HORIZONTAL  ***•«  ************ 

1610  VIEW  ( 60 , 75 )  --  ( 2 1 0 ,  15) 

1  620  W I MDOW  C  -MAX  Z  ,  -M A  X  Z  )  -  C  M A  X  7. ,  M A  X  Z  ) 

1 630  NX  1  ~M I NZ :  NX2--MAXZ :  NY  1  --MINZ :  NY2~ MAX  7. 

1640  LINE  C--NX2, 0)  --  CNX2, 0) 

1 650  L I  ME  ( 0 ,  --NY2 )  -  ( 0 ,  NY2  ) 

1660  REM  ******  ******  DII’U  JO  PRIMER  GRAF  ICO  ***  *  * *****  x  *  »h  x  it  ** 
1670  FOR  I  “  Y  TO  N 
1600  NX  C I ) =C ( I ) 

1690  NY(I)=B(I) 

1700  ’ (  PRINT  C(I),  DC  I) 

17 t0  NEXT  I 

I  730  I .  I  ME  C  NX  <  Y  )  -  0 ,  NY  (  Y  )  -0)  ~  ( N X  (  Y  >  •+  O ,  NY  (  Y  >  +0  >  ,  ,  DF 
1730  ’PRINT  "*" 

1740  FOR  P  Y  i  1  TO  N  129 


1 750  LINE  CNX Cl -1)  ,  NY  Cl'-- 1  >  >  -  C  NX  C  I  >  ,  NY  CIV) 

1760  NEXT  I 

1770  REM  ********  PROYECC I  ON  VERTICAL  PARALELA  AL  A  7.  ******* ***** 

1 780  V I EW  C  250 ,  75 >  --  C  400  p  1 5  > 

1790  WINDOW  C-MAXZ,  --MAXZ  ) -- CMAXZ ,  MAXZ) 

1 800  NX  1 -M I N : NX  2=MA X Z : NY 1 -M IN: NY2-MA X  Z 
1810  LINE  C--NX2,  0)  -  CNX2, 0) 

1 020  L I NE  C  0 ,  -NY2  )  --  C  0 ,  NY2  > 

1830  REM  ************  DIBIJJO  SEGUNDO  GRAF  ICO  ************** 

1840  FOR  I  -  Y  TO  N 
1050  NX  C I ) -PH  C I > 

I860  NY  Cl) -AC  I) 

1070  NEXT  I 

1880  LINE  (NX  CYI-  Q,  NYCYI-Q)-  CNX  C  Y)  '-0,  NY  CY) -i-Q)  ,2,  BF 
1090  FOR  I--Y+ 1  TO  N 

1900  LINE  CNX  C I  - 1  !>  ,  NYC  ID  )-  CNX  Cl)  ,  NY  C  1  )  .) 

1910  NEXT  I 

1920  REM  ************  PROYECC I ON  VERTICAL  PERPENDICULAR  AL  A2  ******** 
1 930  V I EW  C 430 , 75 ) - C 570 ,15) 

1940  WINDOW  C-MAXZ, -MAXZ )- CMAXZ , MAXZ) 

1 950  NX  1  "-M I NZ  :  NX2=MAXZ :  NY  1  =M I NZ :  NY  2- MAX  Z 
1960  LIME  C--NX2, 0)  -CNX2, 0) 

1 970  L INF.  C 0 ,  -NY 2 )  -  C 0 ,  NY2 ) 

1980  REM  ************  DIEtUJO  TERCER  GRAF  ICO  ************** 

1990  FOR  I  -■=  Y  TO  N 
2000  NX  C I ) =PV C I ) 

201.0  NY  C I )  -  A  C  I  > 

2020  NEXT  I 

2030  L I NE  C  N X C  Y ) -Q , N Y  C  Y ) - Q ) - C NX (Y J +Q , NY ( Y > +0 > , 2 , BF 
2040  FOR  I  -Y  i-l  TO  N 

2050  LINE  C NX  C I - 1 ) , NY C I  - 1 > ) - C NX  C I ) , NY C I ) > 

2060  NEXT  I 
2070  LOCATE  21,2 
2080  GOTO  690 
2090  END 


« 
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AIMEXiJ  N°  2 


*  * 

*  RRUcRAMA  "  !•  I  F  "  I  I  l  "I"  i  n  Ik  1 1  P  1 PAbi 

lit'SfE  PPOGRAMA  C  AL.CUl._A  LA  TRANSf 'ORNADA  DC  KMK  i  DA  Ul  LUL/I-m  I: « 

*  CORR1GE  I  AS  AMPLITUDES  Tptos  Uiqit  isadc«  por  minni  >  muJ.  Li*  -n 
ti  (cr.inni  pm-  nn  i  fii  mi-;;  iii.i  i  i  < «  m  nvr,  i  iiut.i  ::  i-i  ai  i.m  i« 
hum  SI.JMV.I  ZADU  1)1:.  HJIlmS  I  I  1M. I  HAN DO  ( il  GU  1/1:  ALIAS  I  A  I  ,'l  H  Ml  I  A  * 


*  Y  CGRRECC I  ON  DE  LA  RESPUES  IA  INSTRUMENTAL  * 

-ft-  !<-  i!  ft**-*  •*  *  *  -.1  H-  •*  *  *  -ft  *  *  -H-Kil  II-  -ft-  ft-  II-  -ft  -ft-  •*  -ft  ft  K  -ft  -ft  -ft-  ft-  -ft  -ft  -k  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft  ft-  -ft  ft  ft  ft  -ft-  ft  ft  ft-  -ft  ft 

ft  EL  PROGRAM  A  I MPR 1  ME  EM  IHPRESOkA  IDS  VALURfcS  LALUJLnniJb  -ft 
ft  C  FRECUEIMC  f  A ,  PER  I  ODD ,  AMPL 1  TUD ,  POTE.NC 1 A ,  I  ASP  C  r  ad  >  .)  .  + 

EL  I  •ROLiRAliA  TAMBIEN  GRAB  A  ElM  ARCH  I.  VO  Sc.CUfc.NC  1 AL  LA  IRE  -  > 
ft  CUEMCIA  Y  AMPL  I  I IJD  PARA  UN  ROSTER  1  UR  l-  LU  I'fcU  POR  HI-.O  1 1 1  ft 
ft  DEL  PROGRAM A  EbPEC.  -ft 

I* 


ft  >  ft  -k  ft  ft  -ft  -ft  -I*  -ft:--ft  K-ft -ft-  x  ft  ft  ft-  hhhhi  ft  -ft  ft  -r;  k  r.  n  **■■**■■*■#■  *  i  h  ■*  +•  >.  -<■  -*  >.  <  ■  r  -ft.  <■  --  r:  ■<•  r.  ^  «.  r.  ■> 

D I  MENS  I  ON  A  C  450©  )  ,  FT  <.  2200  )  .  PH  .1.  •.  22UW .) 

D I  MENS  I  ON  LP  <  4500 )  .  I  P  .1.  ( 600  ) 

COMMON  A,  N,  DT  ,1 
CHARACTER*64  B , F  F I 
CHARACTER  *64  f-FO 
OPEN ( 5 , F I LE” ' CON " > 

OPEN  C  6, 1-  ILE-’CON'  ) 

01  'EN  C 1  ,  E I  LE:-~ '  PEN'  > 

LEO  5 

I MP--6 
4  MM  ~(3 

NR I  TEC  IMP,  700) 

700  FORMAT  C'  T.  DENT  1 F I CAC I  UN  DEL  EVENT  O'?'  > 

READ  CLEG, *  (ACM) * )  B 

WR I  I E  <  1 ,  *  >  B 
WR I TE C I  MR , 70 1 > 

701  FORMAT C'  I NTERVALO  DE  T TEMPO, FACTOR?  ’ ) 

READ  CLEG,  *>  DT, FACT 

WR  1  f  t  (.  1 MP ,  703  > 

703  FORMAT ('  BAJA  Y  ALT A  FREuUENC I  A?'  ) 

READ C LEG, *)  FL.FH  4 

FREGl-FL 

FREQ2TH 

KE=  1 

KG  ~  1 

l-JI-TT  EC  IMP,  707) 

707  FORMAT  C’  N  OMBRE  DE  ARCH  I  VO  DE  DAT  US?1’) 

READ C -ft  ,  ’  ( A) ’  )  EF  T 
OPEN  <7,  FILETFI  ) 

RE ADC  7,  1001)  N 
1001  FORMAT  CIO) 


131 


READ  <7,  1000!) 
1000  FORMAT (2014) 


c  lp  c  i  :> 


i  =  l ,  n  :> 


C*****K3=-1  u  0  P A R A  U N A  ClJI'l I  ••  U I  I E hi  1 E  Y  K 3  =  i  p ,• , i ■ •  ■  t , 

wri  rec  imp,  /uc,:> 

705.  FORMAT*'  NLJITERG  UE  ouiTPONEI'i  I  p.b -  i.U.l  .< 
READ  *  LFO ,  *: ;>  K3 
WRITE*  !  ,  )  T  O 

IF  C K3>  2,2,3 
3  M=N/3 
wO  TO  b 
2  M-N 


5  WR I  TE  Cl,*  >  N 
WRITE (IMP, 938) 

980  FORMAT  ('  PUNTOS:  I  NIC  I AL.  Y  FINAL?  *  > 

RE  ADC  EEC,  *>  N1 ,  N2 
IM=N2-N  1 


J=.l 


DO  955  I  -  I  I 1 , 112 
LP  I  C  J  )  -LP  *  I :» 

C  WR  I  TE  C 1 , 2 .1.  :>  LP  1  <  J  > 

J=J+1 

21  FORMAT <61 4  > 

955  CONTINUE 

WR  I  TE<‘ 1,21 )  C  L.P  ICO.  I  =  1 .  N  > 
C* ****CORRECC 1 ON  POR  EL  FACTOR 


I  PI  .V)  n.u 


DO  99991=1 ,N 
A  C I T -F  LOAT < LP 1 C I > > 

9999  AC  I ) =A  ( I >  *FACT 
PI =3. 14159265 

C  k  *  **  INICIA  LA  REAL  I  ZAC  I  ON  DL  LA  DER1VM  I . UTfiiEK*  •* 

S:=A  C 1  > 

DO  999.1  I  -2 ,  N 
S=S+A C I ) 


9991  CONTINUE 

e>~  s  /  r  i 

DO  9992  1=1, N 
A  C 1 1  -  A  C I  )  -S 

9992  CON  T I  Ml  IE 

C  «•  i*  M  *EJ CUC 1  UN  DEL  OUAV I  Z  ADO  « ■*  *  *  v  t  *.  n 

A  C 1 )  C  A  C 1 )  +  A  C  2  >  >  /  2  ^ 

A  C  2 T  A  C I  :>  +2*A  '•  2  )  +  A  C  3 !>  >  /  4 
Nn=N--2 

DO  9993  1=3, Nn 

A  C  I )  -  C  A  C I  -  2 1  +■  2  *  A  C 1 1  )-i- 3  *  A  (I)  i  2  *  A  C  .!.  -i  O  ■  -  A  ■.  I  +  2  •  ,i  / -j 

9993  CONTINUE 

A  C  N - 1  =  C  A  C  N  -2 1  +2* A  C I  I  -  1  )  +fi  C  N  ,i  >  /4 
A  C  N '  =  < A  C  N- 1 ) -t  A C  N > > / 2 
M--DT  *N 


I.HTI  l  l  I  I- 3 
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WEI TEC  1,12)  DT 

wr  1 1  e  <:  i ,  945  >  c  a  <.  i  :> ,  i  =  i ,  w .» 

94-:,  FORMAT (6E 13. 43 
UR  I  TEC  IMP,  30ICO 

3130  FORMAT <»  NOMBRE  FILE  SAL  1  DA?’  3 
READ'  *,  »  CA.)  ’  >  FFO 
OREM <8,  FIUE-FFQ3 
WRIT  EC8, 301 >  N, Ml , N2, D I 
5(31  FORMA  T  (31 4,  E  10.53 

UP:  J  T  E  C  8 , 900  3  C  A  ( 1 3  ,  T  |  ,  H  > 

•  'OH  FORMA  I  (  1  X  ,  OF  15. 4  3 
■:  CALCULO  Dt.  LA  TPANSF  Ul  MADA 

UP  I  TE (1,10 ) 

CALI.  F  F  I  ( I-  RF  A 1 ,  FREQS ,  I  .  1 ,  F.2 ,  f  T ,  F‘l  11,1  E ,  DF  3 
WRITE( IMP, 863 

86  FORMA  I  '  1.11,  *  E5PEC7R0  DE  OEMS .1  DAD  DE  EMER'.i.lVC' // .) 

7  7  CONriMUE 

T.F  (KG  3  191 , 1.92,  191 

1 91  WRI IE ( IMP, 793 

79  FORMA  I  C  LOGARITHM  SCALE’ // 3 
C  CALL  ESPE C  K 1 , K2 , DF . FT  3 

192  CON  I  l  Ml  IE 

10  FORMAT  ( ’  F  1  LIE ADO  V  AM  Al.  ISIS  DE.  FOURIER  USAMDU  SUBEUIIMA 
1  COOL  ADAP  FACTION  1973’  3 
1 2  FORMAT  C  ’  'I  I  ME  STEP  ’  ,  I  8 . 3/  /  /  3 
WR  I  TE  ( 1.  ,  1. 8  3 

1  fi  EOF  MAI  <  /  /  I  OH  F  1 N  DEL.  AI  L  M.  .1  S  .1  S  3 
IF (MM. EQ. 1 3  GO  TO  4 
S  I  OP 
END 

SUBROUT  I  ME  IFF  ( FREQ  1 ,  FREQ  2 ,  K 1 , 1  2 ,  F  T  ,  PFI 1  ,  Ft ,  bl  3 
D I MENS  I  ON  X ( 4500  3 , I D  C  6  3 , Y ( 1 8 3 
1.)  I  MENS  I  DIM  FT  ( 2200  3  ,  PH  I  C  2200  3 
COMMON  X , N , DT 

C  1 1  IIS  PROGRAM  COMPUTES  THE  FOURIER  .1  Ml  ECiKWL  OF  n  F  UNC  I  .1  UN  ACCURbl 

C  THE  FORMULA  F#X (T>  ZB INTEGRAL  OF  IF  X  Cl  .>  HLXP'.  l*2*FT*f*  I  3  2D  I  llblNu  I  HE 
C  FA 31  FOURIER  TRANSFORM  ALGORYTHM  . 

C  IMF  FOLLOWING  VARIABLES  MUST  BE  GIVEN 
i?  M  ■  NUMBER  OF  POINTS  OF  THE  ARRAY 

C  X  -  X  ARRAY,  WHERE  X  IS  REAL. 

C  I.)T  -  TIME  INTERVAL  * 

C  FREQ1*  FIRST  FREQUENCY 
C  FREQ2*  LAST  FREQUENCY 

C  THE  FOLLOWING  PARAMETERS  ARE  COMPUTED 
C  i:.1.  =  INDEX  CORRESPONDING  TO  FREQ  1 

C  K2  ~  INDEX  CORRESPONDING  TO  I  REQ2C  DF  ~  DELTA  F  REOUENCv' 

C  X'  J  >  ~  THE  ARRAY  OF  MODULI  FROM  K.l  TO  K2„  INCREMENTED  BY  STEPS  OF  2 
C  IE.  J"3 , 5 , / , . 
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C  X<jj:>  =  THE  ARRAY  OF  PHASES  FROM  K1  to  K2,  ll'h  Pfcl'IENIED  BY  S  TEES  Oh 
C  IE.  JJ~J 

C  It.  JJ=4.6,8 . 

C  CHECKS  NUMBER  OF  POINTS  TO  PE  2  10  T  Iil.  I!. 

C  MOTE  BENE  PROGRAM  DESTROYS  THE  DA  I A  AMD  RE! 'I. ACES  .El  WITH  MODULI  >S  •  1  l 
C  AND  PH ASE ( I E . X  C  J  J ) ) 

I  MP-e, 

LEC=b 

P I  -=3 .  1 4 1 5S26535S90 

wr  i  t  e  «:  i ,  *  :>  n 

NMPTS=N 
AN-NMP IS/2 
N0P=8 

L I  =4 

3  L  I  =L OK 
1  =  1  +  1 

IECN.  L..E.  8)  GO  TO  66 
LL--N/L I 

IFCLL.NE. 1)G0  TO  3 

4  N0P=2**I 

IF  (NOP.EQ.N)  GO  TO  5 
1  =  1  +  1 
N0P=2*N0P 
66  NP=NOP~N 
N 1 =N+ 1 

DO  6  K=N 1 , MOP 
6  XCK)=0. 

WRITE'..  I. ,  77)  NP 

77  FORMAT-'  EL  NUHERU  L)E  DATOS  NO  ES  LIN  ITULiiPLO  DE  2  A  LA  N '  ,  /  , 
l'SE  HAN  ANADIDO'  .  17,  '  CEROS  AL  PINAL  DE  LA  HER  IE  PARA  Curi. ’I.ETAk  r  > 
■31.0  FORMAT  f.  3E  13.4) 

5  NE=I-1 
NN= MOP/2 
FNN=NN 
FNOPMTUP 

DF=1.  /<:FN0P*DT)  4 

FNYQ--0. 5/DT 

WR I TE ( 1,13)  NOP , D  T , FNYO , DF 
13  FORMAT  ('  POINTS  NUMBER  =',16,//, 

2’  INCREMENTO  DE  T TEMPO  =',F6.3// 

3'  FRECUENCIA  DE  NIQUIST  =’  F6.3// 

4'  JN1ERVAL0  DL  FRECUENCIA  ='F9.4//) 

CALL  I  '!  COOL  <  ME  ,  '/. ,  D  1  ) 

K 1  --FREO 1  /DF+ 1 
WR I T  E C I NP , * ) K 1 
FI  1  =: ! :‘.  1 

FR='T-  Kl-2.  >*DF 
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IF  CFREQ2.  LT.  FNYO:>  GO  TO  55 
FREQ2-F NYQ 
wr  i  r  i£  < l ,  08  ;> 

QB  FORMAT ( ’  FH  >  NYSQIJIST,  80LU  SALOP AN  FIH  i > 
55  CON  f I NUE 

K2=FREG2/UF H 
7~0.  0 

DC)  7  1  "I!  1,12 

J  J=2* I 

A--3QR1  <X(J>*X(J>+X(JJ:»  *X<JJ)  > 

A~A/AN 

IF  CX«:j:>  >  97 ,96,97 
96  B~P I / 2 . 

GO  TO  98 


HR  NY  SOLUS  I 


9/  CONTINUE 

St:J0  FORMAT  (  ’  TYPE  ITAGN  IF  T  TAT  I  ON  I-  (IF:  IITIS  SEJ  SMOGRaPH  lUl'IL  Y  ON' 
8- AT AN 2  C X C J J > , X C  J > > 

98  Cf  I  N’T  I  NUE 

IFCB.LT.e.  )  B~B+2 .  *P I 

ft  <:  i  :>  ~a 
phi  <  i  :>=b 


7  CONTINUE 

IF(KE)  32,84,83 
03  CALL  RESP ( FT , K 1 , K2 , OF  > 

GO  TO  84 
02  CONTINUE 

WRITE (IMP, 500) 

READ (EEC,  *)  PEAK 

21  FORMAT  (F"  10. 4) 

WRITE  U,  22)  PEAK 

22  FORMAT  < »  VALOR  CORREGIUO  FUR  IMS  TRUITENT  O.  AMI-'L.  MAX"  ’,18. 
DO  27  I=K1,K2 

FL=(I-1 . )*DF 

ftci  :>  ~f  r  a )  /  se  i  sma  c  fe  ,  peak  .> 


27  CONTINUE 
84  CONTINUE 
WRITE!.  1,  10) 

1 0  FORMAT  <//»  I  TIDE  X  ’  4  X  ’  FREQ  ’  4  X  ’  PER  l  UD  ’ 
1 » p H A 3 E ( R AD)  P H  C U M ' / /  > 

R=0 . 

M 1 "K2-K 1 + 1 
WRIT  EC 8, 901)  Ml 
901  FORMAT C 14) 

DO  &  I  ~K  .1  #  K2 
PS  -  F  r  ( I  >  *F  I  ('  I ) 


5X»  AMPE’ S) 


P0TENC1 A’ 


4 


Z-2  +  RS 
PK-PH  HI) 
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IFCI.EQ.K1)  GO  TO  1(31 
IF(PK-PKB) 100,101,101 
10(3  R=R+1. 

101  PHI  (I  !>  -PH  I  Cl)/  C  2 .  *P  I  :>  +R 
F'R-FF'+DF 
PKB=F‘K 

IF  CFR)  90,91,  S'/l 

91  FT— 9999. 9999 
GO  10  92 

90  TT-l./FR 

92  WR II  E  C  1 ,  1 1  )  I ,  FR ,  TT ,  F T  C  I  )  ,  PS ,  PL ,  PI  1 1  «. 1  .» 

1 1  FORMAT  Cl  5 ,  2F 1 0 .  4 , 2E 1 2 . 5 .  2F 1 0 .  4  .> 

UR  I TE C  & , 900  )  FR , FT C I ) 

900  FORMAT C FI 4 . 4 , E 1 7 ; 5 ) 

S  CONTINUE 
9  CONTINUE 

Z  Z-SQRT C  Z* . 0004327 ) 

WRITE C 1 , 40>  Z,Z.Z 

40  FORMAT  C f  LA  PUMA  DE  LOS  Cl.li  ibRPiJHJ'  >  UL  t.UH  MOut t:Sr  ,  I.  I  I ,  /  , 
1 '  EL  LIMITE  DE  CONE  I  AN/ A  DEL  95  POP  C1ENI0  ES’  .L12.G//J 
19  CONTINUE 
RETURN 
END 

SUBROUT  I NE  RECOOL  C  N ,  X  ,  D  T '.) 

DIMENSION  X  C  1  >  ,  6 C 2 .) 

COMPLEX  X , A, B, W, CON  IG 
EQU I VALENCE  C  6 , W > 

D2=0.5*DT 

L=2*-xN 

WRITE  Cl,*)  COMPLEX  X 
CALL  COOLCN, X, -1 . > 

FL=L 

ARG=3. 14 159265 /FL 
G C 1 ) “COS C  ARG  > 

GC2)-SINCARG) 

WR IT  E  C 1 , * >  X  C 1 ) 

B—CONJG  C  X  C 1  )  :> 

A-  X C 1 ) 

WRITE  Cl,*)  B , A 
X  C 1 )  =  C  A+B+ C  0 . , 1 . ) * C  B-A > ) *D2 
X  C L+ 1 ) - C A+B- C 0 . , 1 . ) *  C B-  A ) ) *D2 

LL-L/2+ 1  4 

DO  10  I "2, LL 
J-L  —  I  t-2 
B-CONJG C  X  C I ) > 

A--B+X  C  J  .) 

B~ C  X  C  J ) -B ) *W** C I  —  1 ) 

X  C  J )  =  C  A + C  0 . , 1 . ) *B ) *D2 


-i 
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10  x  <  i  > - 1  cdnjg <: a :>  +•<: 0 . ,  i .  :> *con j g <  b  >  *  *i>,- 

RET URN 
END 

RJNCT  I  Dl'l  SE  I  SNA  <  FE ,  FEAK  .> 

ESI  A  FUNGI  ON  CALCULA  LA  RESPULb  i  IWb  IKUMLNi  Ai. 

PARA  l ..AS  ESTAC30NES  DE  LA  RED  WCJ3S 
WE=6. 2331853*FE 
IMP =6 

1 2  FORMAT C 2X , » FE= ' , F 1 0 . 5 . ' SL I SMA= * , F 1 0 . b > 

I F  < PEAK . EO . 375 .  )  GO  T O  1 
IF '..PEAK.  EO.  750.  >  GO  10  2 
I F  <  PEAK .  ED .  1 500 .  '.>  GO  TO  3 
I F  C PEAK . ED . 3000 . >  GO  I  0  4 
I F  C PEAK . ED . 6000 .  >  GO  TO  5 
t  FMAG=298. 44574 
SIGMA-0.003 
GO  1 0  6 

2  F  M  A  G  -  5  9  6 .  4  /  3  3  0 
SI GHA ■ 0.013 

GO  TO  6 

3  I  'MAG-- 1 1 89 .  99848 
SI  GMA--0 . 04  / 

GO  TO  6 

4  FI1A6-2349.  20/90 
S I GMA~0 . 204 

GO  TO  6 

5  FNAG--4236. 86538 
S I  GI4A~0 . 305 

6  ZETA~0 . 93 
2 ETA 1  ~  1 . 0 
WN~0. 4188/9 

WN 1=0. 062831853 

AR~  C  WE  *WE-WN  *WN  !>  *  <  WE*  WE-WN 1  *WN  1  > . 4 . 0  *2E  T  A  *  ZE  I  A 1  *WN*WN  1  *  C  l .  0  -Vi  I GMA 

1 ) * WE* WE 

A I -2.  0*WE*  (.  ZETA1  *  WN1  *  {  WN*WN--WE*WE )  +ZE  I  A  *WN*  (.  WN  1  *WN  1  -WE  *  WE  .> .» 

SE 1 3I4A~FMAG*WE*WE*  WE /SORT  <  AR*AR+  A I  *A  I ) 

WR I T E C I MP , 1 2  >  ( FE , SE I SMA > 

RETURN 

END 

3LJBR0UT I NE  COOL  C  N ,  X  ,  SW ) 

D I  MENS  I  ON  X  (.  1  :>  , I NT  <  1 6  >  ,  G  <  2 !> 

COMPLEX  X , Q, W, HOLD 
EQUIVALENCE C 6, W i 

LX-2**N  * 

FLX-'LX 

IL-LX 

PI. 2=6. 203105306 
i-  L  X  P 1 2  •=  S  W  *  P 1 2  /  F  L  X 
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DO  10  1  =  1  ,N 
IL=IL/2 
10  I  NT  C I  :>  =  I L 
NBLOKK 1 

DO  40  MATE I  X " 1 , N  . 
NBLOCK=NBLOKK 
NBLOKK=NBLOKK+NBLOKK 
LbLrjCK =L  X  /  NBLGCK 
LBH  ALR  =  L  BLOCK  /  2 
NW=0 

DO  40  I BLOCK- 1 , NBLOCK 
LST ART  =L  BLOCK*  (  1  BLOCK- 1 
FNW=NW 

AR6=FNW*FLXF‘  1 2 
Q< 1 ) “COS  C  ARG> 

6  f.  2 >  ~S  I N 1  ARG 
DO  20  1  =  1,  LBHAL.F 
J= I +LST ART 
K= J+LBHALF 

c  wr  n  e  <  i ,  * >  <:  x  c  i  :> ,  i = 1 ,  lbhalf  > 

G=X(K>*W 
X ( K )  =  X < J)-Q 
X  < J)=X (J  >+Q 
20  CONTINUE 
DO  32  1=2, N 
IK- 1 

I F  (.  I  NT  <  I  >  -NW :» 32 , 32 ,  40 
32  NW=NW-INTCI> 

40  nw=nw-hnt<:ik:> 

NW=0 

DO  80  K=- 1 ,  LX 

nwi=nw+ 1 

IF(NWl-K)  55,50,60 
60  HOL.D=X  CNW1  > 

x  c  nw  i  :>=x  oo 

X  ( K  :>  =HULD 
55  CONTINUE 

DO  70  1=1, N 
IK  =  I 

I F  Cl NT < I > -NW  >70,70,80 
70  NW=NW- 1 NT  C I > 

80  NW--NW+INT  C I K  > 

RE  TURN 
END 

FUNCTION  CON  J  Li  C  L  > 

L)  I  MENS  I  ON  2  C  2  > 

COMPLEX  CON  JG 

C0WJG=2  C  1  >  -  c  0 .  ,  1 ,  :>  *7.  ( 2  :> 
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EL  TURN 
END 

St  IPROUT  1  NE  ELSE  <  I  f  ,  F,  l ,  1,2 ,  OF  ) 
i:  HLIL  FI  PLICA  LOS  YALOEEb  DLL.  I  SI  t  i  1 1:0  lit  ni  in 

C  DR  LA  CURVA  UL  RESPULGTA  DLL  lllb  I  Pr.l  il  l'll  H  (i: 

C  -LOS  7AL0EES  FUEPri  DE  I  A  L  I  iRVi .  BE  U  F  1 1.  I  .  A.  I 

C  PUN  10.  LOB  7AL.UEES  DE  L.  A  L  UkVA  l;L  Kt- OPUES  I , 

C  RESPECTO  A  PERIODGS  Y  SE  ALMACLNAN  EN  X  U  > 

C  Y  AMPLITUDES  EN  X  ( 2  *  ,T- 1  )  . 

C  NPA— NUMERO  DE  PUNT  OS  EN  LA  CURVA 

C  VMT  S— VALOR  DE  l.A  CURVA  PARA  EL  VALOR 

C  NORMAL I Z ADA 

C  X  C  J  :>  -CURVA  -  NORMAL  1  Z  ADA 

D I  MENS  I  ON  X  ( 4500  >  ,  E  I'  ( 2200  > 

COMMON  X ,  N , DT 
CHARACTER*64  RAO 
OPEN ( 5 , F I L  E= f  CON ” ) 

OPEN C 6, FILE-' CON'  > 

□PEN  <  1 ,  F I  l.E~  ”  L.P  T  .1 "  .) 

READ  C 5 , 28  >  NPA , VM IS 
N2PA=2*NPA 
N I  PA=N2PA-  .1 
UR  r.  IE'  [MP,  /in) 

/as  FOF'LMA  r »'  1H,  *  NOM£«RE  DE  LA  lUEVh  Dt  kk.Bi -ULS  i  A 
READ'*, » LA) ” >  RAG 
OPEN  <  9 ,  F  I.  LE=R  A'3 ) 

REAL'i  9,  *>  C  X  (  I  )  ,  1  =  1 ,  N2PA> 

DO  500  I  "•  1  ,  N2PA 
500  X  C  I )  -  AL.OG  1 0  (.  X  (  I )  ) 

DO  31  JO-Kl ,  l<2 
F  JO=JO 

FR=<FJ0-1 . >*DF 
IF  (FR) 20, 2 1,20 

21  rr= loooo.w 

1  iO  I  0  22 

20  rr=i./FR 

22  CONTINUE 

r  I-- a  log  it)  r  r;> 

DO  I 2  J  = I , Ml PA, 2 
L---J  + 1 

IF  O  I  —  X  (L.)  »  1.0,  11,12 
10  IF'  (  I  I  -  X  >2')  >  13,  14,  14 

1 3  77  X (  1  ) 

GO  r  0  15 

14  DIFT“X<L)-X<L-2) 

D I  F'7~-X  <J)-X  <  J  —2  ) 

DIF'-T  1  -X  0.-2) 

W=CDIFV/DIF7  )*DIF+X  <  J-2» 


.j  lliL'tb  F'uE  I  n  j  A  XhJ-'.'.i  1 
:<A  1  N  I  LRPULAC  L  Al  l  L.  I  Ml::  l  iL. 
EL..  RRU-IF..P  Y  NI..  lira  I 
T  EB  I'AN  LEI  l)U'B  LUN 
--F  ER  101)0S  EN  X  '  2*  J  .» 


RESPECTO  nl.  t!UL  ESI  A 


'  .) 


I 
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Cr*  CT* 


GO  TO  15 

11  W-X<J> 

15  W=lu-*VvV 

f  r  (  jo  :>  ~  r  r  <:  jo  '» /  <:  vv*vi  ri  s  .> 

GO  TO  31 

12  CONTINUE 
VV=10**X  C I  I 1  PA  ) 

r  r  c  jo  :>  -f  r «:  jo  :>  /  <:  vv  *  vmts  > 

31  CONTINUE 

WRI  TE  ( 1 , 25". > 

25  FORMAT  < » LOS  VALQRES  ESTAM  CORREGIDOS  PUR  LA  CURVA  uEl_  li  IS  VRUi-lh  l  l  I  u" 
l//> 

DO  26  J=1,N2PA,2. 

L=J+1 

x  <  j  :>  =  x  c  j  :>  +al  oo  1 0  <:  vm  s  :> 
x  c  j:i  =i0**x  <  j:> 
x  (l:>  =  10*!'-x  <:l:> 

CON T I NUE 

WRITE  <6,  27)  XCL)  ,  X(J) 

2 /  FORMAT <5X,F10.5, A X.F1 2 . 2 > 

28  FORMAT C 1 3 , FA  -  2  ."> 

RETURN 

END 


4 


> 
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ANEXO  N°  3 


8  • 

10  REM  x  -  — . .  PEOGEAMA  LSPLC  - - — . - - -  * 

20  REM  *  ESTE  PEOGEAMA  DIBUJA  LA  ON DA  Y  SU  ESPECTRO  DE  * 
30  REM  *  AMPLITUDES  CALCUL.ADQ  POE  EL.  F  AST-FOUR  I  EE.  * 

40  PPM  *  -  El...  PEOGEAMA  HA  SI  DO  ESCRITO  EN  GWBASIC  PARA  * 

*  50  REM  *  I..A  COMPUTADORA  ACER-  710  DEI  OPS.  SAM  CAL  I XI 0 . -* 

G0  REM  *  OCT,  1389  E.  MINAYA.E  * 

70  REM  *  - - -  * 

f  30  REM  *  PARA  OBTENER  LA  IMPELS I ON  DEL  GRAF I CO  SE  * 
90  REM  *  SIGUEN  LOS  SIGUIENTES  PASOSs  * 

100  REM  *  C  >CD  DOS  * 

110  REM  *  OGRAPHIC5  * 

130  REM  *  _  Oft:  * 

130  EEM  x  A > GWBASIC  * 

140  EEM  x  OK  * 

150  EEM  *  LOAD" ESPLC”  * 

1G0  EEM  x  UNA  VEZ  QBTENIDO  EL  GRAF  ICO  SE  PRES IONA  LA  * 
1.70  EEM  x  TECLA  IMF’EIME  PAN!  ALLA  * 


190  OPTION  BASE  1 

200  DIM  A(  1024)  ,  TEC  1024!) 

210  DIM  MX  <  750 >  ,  NY  C  750 )  ,  EE  <  750 !>  ,  FT  C 750  > 

220  ’  INPUT  "  I  DENT  I E I CACI  ON  DEL  EVENTO  "  ,  NAM'f 

230  INPUT  "ARCH I VO  " , I F I  LET 

240  INPUT  "POSICION  GRAF I CO?" , XX 

250  OLS 

2G0  SCREEN  2 

270  OPEN  "I",  1,  IF  I  LET. 

200  INPUT#  1 , N , DT 
290  FOR  1=1  TO  N 
300  INPUT#  1 , A C I > 

31.0  NEXT  I 

320  TPM-DT*f.N-l  > 

330  IT  XX- 1.  THEN  GOTO  3G0 

340  IE  X X  ~ 2  THEN  GOTO  370 

350  IE  XX =3  THEN  GOTO  300 

3G0  T  -0:  NX2-430:  NY  1 --80:  GOTO  390 

370  T-0: NX 2=250: NY  1=80: GOTO  390 

300  T-0 : NX2=G0 s NY  1 -80 : GOTO  390 

390  PEM 

400  PEM  *********  DIBUJA  F-JES  DEL  ESPECTRO  ****** 
■110  PEM 

420  l  IMF  <NX2,  NY1  :>--<NX2  i  B0,  NY1) 

430  LINE  (NX 2, MY  1 )  - CNX2, NY 1 -601 
440  FOR  1=1  TO  4 
450  NX2-NX2+20 

4  80  LI  ME  (  N  X  2 ,  NY  1  - -2  >  -  <  M  X  2 ,  MY  1+29 
470  NEXT  I 

400  IE  XX- 1  THEN  NX 2-430 
490  IP  XX -2  THEN  NX2-250 
500  IP  XX -3  THEN  NX 2-80 
510  POP  1-1  TO  4 
520  NY 1 -NY 1-15 

530  I  I  ME  <  N  X  2  4 ,  NY  1 .»  -  ( N  X  2 + 4 ,  NY  1  !> 

540  NEXT  I 
550  NY  t  —00 
580  PPM 
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570  REM-*****  DIBUJO  DEL  ESPECTRO****** 

580  REM 

590  INPUT#  1 ,  M 

600  FOR  1=1  TO  M 

610  INPUT#  1 , FR  < I ) , FT  ( I ) 

620  NEXT  I 

630  AMAX=FT  (.  1  >  :  AMIN=AMAX 

640  FOR  1=1  TO  M  ' 

650  T=FT  < I > 

660  IF  T>=AMIN  THEN  GOTO  630 

670  AMIN=T  « 

680  IF  TOAMAX  THEN  GOTO  700 
690  AMAX=T 
700  NEXT  I 


710 

IF 

AMAX  O.  0001 

THEN 

T4=— 7 : 

GOTO 

760 

720 

IF 

AMAXO.  001 

THEN 

T4=- 6: 

GOTO 

760 

730 

IF 

AMAXO.  01 

THEN 

T  4=~5: 

GOTO 

760 

740 

IF 

AMAXO  .  1 

THEN 

T4=-4: 

GOTO 

760 

750  IF  AMAXO  1  THEN  T4=~3:G01Q  760 

760  IF  XX= 1  THEN  NX2=430 

770  IF  XX=2  THEN  NX2-250 

730  IF  XX =3  THEN  NX2=60 

790  DF  - 45: DFT= 1 5 : NY  1 =80 

300  FOR  1=1  TO  M 

8 1 0  NX  < 1 1 =NX2+FR < I ) *DF 

820  NY  ( I  )  =NY  1  --  <  < LOG  ( FT  <  I  >  )  /LOG  C 1 0 )  )  -T4+ 1  !  >  *DF  T 

830  NEXT  I 

340  FOR  1=2  TO  M 

850  LI  NE  (  N  X  ( I  - 1  :> ,  NY  <1-1 )  -  C  NX  Cl),  NY  <  I  )  ) 

360  NEXT  T 
070  KEY  OFF 
380  END 


« 


v 
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Translation  of  Thesis  Abstract 


y  National  University  of  San  Agustin  de  Arequipa  (Peru),  Department  of 

Geology  and  Geophysics 

Rq  Waves  fmm_Pemyian  seismic  events  recorded  at  the  La  Paz. 

*  Bolivia  (LPB1  station 

Thesis  (undergraduate),  Celedonio  Mora  Tito 

Arequipa,  Peru,  1990 

Rg  waves  from  121  earthquakes  located  in  Peru  or  in  the 
offshore  region  adjacent  to  Peru,  and  recorded  at  La  Paz,  are 
examined  in  the  context  of  similar  studies  completed  in  other  parts 
of  the  world.  General  geological  characteristics  of  the  area  are 
reviewed,  as  are  the  principal  geophysical  methods  previously 
applied  to  the  source  region  and  to  the  Rg  propagation  paths;  in 
particular  velocity  profiles  between  pairs  of  stations  determined  by 
James  (1971)  are  considered  in  this  review. 

Earthquake  sources  in  this  study  have  been  separated  into 
three  separate  zones  based  on  seismogram  characteristics:  South, 
Center  and  North.  Velocities  do  not  differ  much  between  zones  in 
the  2.8  to  3.4  km/sec  range.  The  first  zone  has  large  Rg  amplitudes, 
without  apreciable  variation  between  different  geological  provinces 
except  for  reduced  amplitudes  from  offshore  epicenters. 

The  second  zone  exhibits  a  smooth  amplitude  gradation  from 
larger  amplitudes  in  the  eastern  part  through  diminishing 
amplitudes  in  the  more  western  geological  provinces  to  negligible 
amplitudes  from  offshore  earthquakes.  In  the  third  zone  the 
amplitudes  are  smaller,  with  attenuation  attributed  to  the 
propagation  path,  and  with  a  gradation  similar  to  that  observed  in 
Central  Peru  (Western  and  Eastern  Cordilleras). 

Computer  analysis  of  particle  motion  is  applied  to  hand 
digitized  records  for  a  variety  of  propagation  paths,  incorporating  S 
wave  decomposition  (this  being  the  factor  that  can  have  the 
greatest  complicating  effect  on  the  analysis).  Characteristics  of 
Rayleigh  surface  waves  are  prevalent  in  the  particle  motions. 

Fourier  spectra  from  FFT  analysis  display  a  cutoff  frequency 
of  0.8  to  1.05  Hz,  above  which  amplitudes  diminish  rapidly.  Spectral 
amplitude  maxima  ranging  from  0.2  x  10*5  mm  to  1.0  x  1 0'5  mm 
occur  at  frequencies  of  0.7  to  0.9  Hz.  Intercalated  minima  do  not 
have  great  significance.  The  Rg  waves  propagate  through  the  crust, 
probably  guided  through  a  layer  the  top  of  which  is  located  5  or  6  km 
below  the  surface  (granitic  layer). 
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